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' Abstract: 

We present a complete description of the spectrum of compound states of reggeized gluons in 
QCD in multi-colour limit. The analysis is based on the identification of these states as ground 
states of noncompact Heisenberg SL{2, C) spin magnet. A unique feature of the magnet, leading 
to many unusual properties of its spectrum, is that the quantum space is infinite-dimensional 
and conventional methods, like the Algebraic Bethe Ansatz, are not applicable. Our solution 
relies on the method of the Baxter Q— operator. Solving the Baxter equations, we obtained the 
explicit expressions for the eigenvalues of the Q— operator. They allowed us to establish the 
quantization conditions for the integrals of motion and, finally, reconstruct the spectrum of the 
model. We found that intercept of the states built from even (odd) number of reggeized gluons, 
A, is bigger (smaller) than one and it decreases (increases) with N approaching the same unit 
value for infinitely large N . 
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1. Introduction 



It has been recently realized that QCD possesses a hidden symmetry at high-energy |T|, 0. 
This symmetry is not seen at the level of classical QCD Lagrangian and it manifests itself 
through remarkable integrability properties of the Schrodinger equation for the partial waves of 
the scattering amplitudes in perturbative QCD in the so-called generalized leading logarithmic 
approximation (GLLA) [0, It turns out that in the multi-colour limit this equation coincides 
with the Schrodinger equation for two-dimensional quantum-mechanical completely integrable 
model, which was dubbed in p, ^ as a noncompact Heisenberg spin magnet. 

The asymptotics of the scattering amplitudes ^(s, t) at high energy, s ^ — t, is governed by 
the contribution of an infinite number of soft gluons exchanged between the scattered particles. 
In the GLLA approximation, the scattering amplitude is given by |P, ^, ^ 

^(S,t) ~ -tsj^i^^sfj-^ 7-^UN{t)^bAi), (1-1) 

^2 (asCTTvlns)!/^ 

where as = c^sNc/tt is the QCD coupling constant and the sum goes over an arbitrary number 
of reggeized gluons exchanged in the t— channel, = 2,3, .... In the GLLA approximation, the 
reggeized gluons interact with each other elastically and form colour-singlet compound states 0]. 
These states can be defined as solutions to the Bartels-Kwiecinski-Praszalowicz (BKP) equa- 
tion I, I 

The effective QCD Hamiltonian Ti^ acts on the colour SU{Nc) charges of A^ reggeized gluons 
and their two-dimensional transverse coordinates, Zk {k = 1, ...,N), which belong to the hyper- 
plane orthogonal to the momenta of two scattered particles. The contribution to the scattering 
amplitude ([Ll| ) of the compound states built from A^ reggeized gluons, "^{zi, Z2, zn), has the 
standard Regge form ~ g-as-Biv/4_ jg dominated at large s by the contribution of the ground state 
En for Eq. (|1.2| ). Obviously, it increases (or decreases) with s if the energy of the ground state 
is negative (or positive). As we will show below, the spectrum of (|1.2|) is gapless (see Eq. 



below) and, as a consequence, one has to retain in ( |1.1| ) the contribution of the excited states 
next to the ground state. This amounts to appearance of the additional factor (asO"Ar In s)~^/^ 
in the r.h.s. of (p^). The residue factors ^a{f)),Af(t) measure the overlap of \E'(zi, ^2, iW) with 
the wave functions of the scattered particles and they depend, in general, on the momentum 
transferred, t, and the colour factor, l/N^- 

Calculation of the spectrum of the compound states (|1.2| ) for arbitrary number of reggeized 
gluons A^ and eventual resummation of their contribution to the scattering amplitude (|1 . 1| ) is a 
longstanding problem in high-energy QCD |l^, |Tl], |T^. At A^ = 2, the solution to ( |1.2|) has 
been found a long time ago - the well-known Balitsky-Fadin-Kuraev-Lipatov (BFKL) Pomeron 
[0]. At A^ = 3 the solution to (|1.2|) - the Odderon state in QCD |[l^, was formulated only a few 
years ago by Janik and Wosiek ||T^ by making use of the remarkable integrability properties of 



the effective QCD Hamiltonian 0, ^. Their solution has been later verified in Refs. [jT5|, |T^. The 



methods employed at A^ = 3 in can not be generalized, however, to higher (A^ > 4) reggeized 



gluon compound states and very little is known about the solutions to (|1.2|) for A^ > 4. Recently, 
a significant progress has been made in solving the Schrodinger equation (|1.2|) for higher A^ in 
the multi-colour limit, ^ oo and = fixed 0, |T6[. The first results of the calculation of 



the ground state energy E^ for higher reggeized gluon compound states in multi-colour QCD 



2 



were reported in a letter format [O. In this paper, we shall provide a detailed account on the 



approach used in [T^] and present new results for the spectrum of the Schrodinger equation ( \L2i) . 
The reader interested in learning more about the latter could skip the first part of the paper and 
go directly to Section 5. 



Our approach to solving the BKP equation ( |1.2[ ) is based on the identification of the effective 
Hamiltonian in the multi-colour QCD, Tijv, as the Hamiltonian of a completely integrable two- 
dimensional noncompact Heisenberg spin magnet [|T], 0. The latter model describes the nearest 
neighbour interaction between spins of particles "living" on a two-dimensional plane of trans- 
verse coordinates z = {x, y). The corresponding spin operators Sk and Sk (with = 1, N) are 
the generators of the unitary principal series representation of the SL{2, C) group specified by a 
pair of complex spins (s, s). They act on the i*— plane as the differential operator^ 

cO 



Zkdz,. + s , Su = -d; 



Si, 



;i.3) 



so that Sk = {SlY + {S^S^ + S^S^)/2 = s{s — 1). The operators Sf'^ are given by similar 
expressions with Zk and s replaced by Zk and s, respectively. Here, the notation was introduced 
for the (anti)holomorphic coordinates on a two-dimensional £— plane, Zk = Xk + iyk and Zk = z^, 
so that d?Zk = dzkdzk/2. By the definition, [S"^, S"^] = for a, 6 = ±,0. For the principal series 
of the SL{2, C), the possible values of the complex spins (s, s) take the form |T8| 



;i-4) 



2 ' 2 

with Us integer and Vs real. For the reggeized gluon compound states, Eq. (|1.2[), the 5^(2, C) 
spins take the values s = and s = 1, or equivalently Us = —1 and z/^ = 0. 

The Hamiltonian of the noncompact SL{2, C) Heisenberg spin magnet is given by [0, Q 



N 



T^N = ^ [H{Jk,k+i) + H{Jk^k+i)] 



:i.5) 



k=l 



where H{J) = ip{l — J) + ipiJ) — 2ip{l) with ip{x) = d\nT{x)/dx, Jk,k+i is the sum of two 

SL{2,C) spins, Jk,k+i{Jk,k+i - 1) = {Sk + Sk+iY with Jn,n+i = Jn,i, and similar for Jk,k+i- 
The model possesses the set of mutually commuting conserved charges qk and qk [k = 2, ...,N). 
Their number is large enough for the Schrodinger equation ( |1.2D to be completely integrable. 
The charges qk are given by the A;— th order differential operators acting on the holomorphic 
coordinates of particles. They have particularly simple form for the SL{2,C) spins s = and 

s = iE§ _ 

:i.6) 



qk 



s=0,s=l 



l<jl<j2<:.<jk<N 



Ik 



with z. 



■jk 



Zj — Zk- The charges qk are given by similar expressions in the sector 



qk 



s=0,S=l 



i<ii<i2<...<ife<Af 



•••^^ji._i'^2ji.^iii2---%-i,ifc^ifcji 5 



[1.7) 



so that qk = ql with respect to the SL{2, C) scalar product 



1^1 



d^Zi...d^ZM \'^{zi, Z2, Zm) 



^That is, the isotopic "spin" space coincides with the coordinate space of N particles. 
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The eigenstates £2, ^Af) have to diagonahze these operators and be normahzable with 

respect to ([L^ ). The corresponding eigenvalues q = {qk,(lk = 5^}, with k = 2,...,N, form the 
complete set of quantum numbers parameterizing the spectrum of the Schrodinger equation (pr2|), 
= Eff[q, q). The eigenproblem for the operators ( |1.6|) and (1/7) leads to a complicated system 
of (A^—1)— differential equations on "^{zi, 22, ^a^), which was previously solved at = 2 ||^ and 



N = 3 |Tj]. For higher A^, instead of dealing with this system, we apply the method developed in 
[^]. It represents an application of the Quantum Inverse Scattering Method to noncompact 
Heisenberg spin magnet model. 

The noncompact Heisenberg spin magnet, Eq. ( |1.5D , can be considered as a generalization of 
the well-known spin— 1/2 Heisenberg spin chain model (as well as its analogs for higher SU{2) 
spins to arbitrary complex spins belonging to noncompact, unitary representations of the 
SL{2, C) group. As we will demonstrate below, the noncompact and compact Heisenberg magnets 
have completely different properties, yet another manifestation of the fact that the quantum 
space of the model is infinite-dimensional in the former case. In particular, in distinction with the 
compact spins, the principal series of the SL{2, C) group does not have the highest weight and, as 
a consequence, the conventional Algebraic Bethe Ansatz is not applicable to diagonalization 
of the Hamiltonian (|1.5| ) . To solve the Schrodinger equation (|1.2| ) for arbitrary number of particles 
A^ we will rely instead on the method of the Baxter Q— operator . 

In this method, the Hamiltonian ( |1.5|) , the integrals of motion ( p..2|) and, in general, all 
transfer matrices of the noncompact Heisenberg spin magnet are expressed in terms of a single 
operator Q{u,u), which acts on the quantum space of the model and depends on a pair of 
spectral parameters, u and u. The explicit form of this operator was found in p. As a result, 
the Schrodinger equation (|1 . 2| ) turns out to be equivalent to the eigenproblem for the Baxter 
Q— operator. It is this problem that we address in the present paper. Namely, we calculate 
the eigenvalues of the Baxter Q— operator, establish the quantization conditions for the integrals 
of motion (|1.6| ) and, finally, obtain a complete description of the spectrum of the Schrodinger 
equation (|1.2[). 

As we will show below, the system (|1.2|) has many features in common with two-dimensional 
conformal field theories (CFT) |2^. Since the Hamiltonian (|1.5| ) is given by the sum of two 
mutually commuting operators acting in the z— and sectors, the dynamics in the two sectors 
is independent on each other. As a consequence, the solutions to the Schrodinger equation 
( |1.2|) have the chiral structure similar to that of correlation functions in the CFT. Namely, the 
eigenstates \E'(ii, ^2, zn) can be factorized into the product of "conformal blocks" depending on 
the (anti)holomorphic coordinates and the conserved charges q. Similar factorization holds for the 
eigenvalues of the Baxter operator Q(m, u). For "^{zi, Z2, z^) to be a single- valued function on 
the two-dimensional z— plane, the conserved charges q have to satisfy the quantization conditions. 
The charges q play the role analogous to that of the conformal weights of primary fields in the 
CFT. As we will show, the spectrum of their quantized values turns out to be very similar to the 



Kac spectrum of the conformal weights in the minimal CFT 



The paper is organized as follows. In Section 2, we summarize the main properties of the 
Baxter Q— operator for noncompact Heisenberg spin magnets. In Section 3, we show that the 
problem of finding the eigenvalues of the operator Q('U, u) can be reduced to solving the A^th 
order Fuchsian differential equation. Its solution leads to the set of consistency conditions which 
can be satisfied only if the integrals of motion q take quantized values. In Section 4, we calculate 
the eigenvalues of the operator Q(m, u) and demonstrate that they are factorized into a product 
of "conformal blocks", which depend separately on the (anti)holomorphic spectral parameters, u 
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and u. Using these expressions, it becomes straightforward to determine the exact spectrum of 
the model for arbitrary number of particles, A^, and complex SL{2, C) spins, s and s. In Section 
5, we present the results of our calculations for the special values of the SL{2, C) spins, s = and 
s = 1. The obtained expressions define the spectrum of the compound states of reggeized gluons 
in multi-colour QCD. Section 6 contains the concluding remarks. The details of the calculations 
are summarized in the Appendices. 

2. Baxter Q— operator 

In this Section we shall describe, following , the general properties of the Baxter Q— operator 
for the noncompact Heisenberg spin magnet. We assume that the number of particles is 
arbitrary and the complex 5*17(2, C) spins s and s are given by (pT^) . The operator Q(m,'u) 
depends on two complex spectral parameters u and u. It acts on the quantum space of the model 
Vn = ® ... ® V^, with V = y being the representation space of the principal series of the 
SL{2,C) group. For Q{u,u) to be a well-defined operator on Vn, the spectral parameters have 
to satisfy the additional condition 

i{u — u) = n (2.1) 

with n being an arbitrary integer. The Baxter Q— operator commutes with the Hamiltonian of 
the model ( |1.5|) and shares the common set of the eigenstates 

Q(m, u) -^p^{q^q]{Zl, Z2, Zn) = Qq,q{u, u) ^p-,{g,g} (^1 , ^2, ^Tv) , (2.2) 

with p being the total two-dimensional momentum of the state and {g, q] denoting the total set 
of the quantum numbers, qk and % with k = 2, N. 

The Baxter operator plays the central role in our analysis as the energy spectrum of the 
model can be expressed in terms of its eigenvalues, Qq^q{u,u). Indeed, there exist the following 
two (equivalent) relations between Qq^q{u, u) and the energy = E^i^q, q) 

EN{q,q)=eN + i^\n[QJu + is,u + is){Qgg{u-is,u-is)y] , (2.3) 

where En = 2NRe [^(2s) + ^/'(2 - 2s) - 2tp{l)] and 

ENiq,q) = -Im^ln u^^ Qq^q{u + i{l - s),u + i{l - s)) (2.4) 



xQ^q^_q{u + i{l - s),u + i{l - s)) 



u=0 



The Hamiltonian ( |1.5| ) is invariant under cyclic permutations of particles [P\I'](zi, ^at-i? zn) = 
'^{z2, zn, Zi) and, as a consequence, its eigenstates possess a definite value of the quasimo- 
mentum defined as 

[P%{,,,|](2;i,...,zjv) =e*'"(^'^") %,,,-}(2;i,...,zjv), O^iq, q) = 2'Kk/N , (2.5) 

with k integer in virtue of = 1. The quasimomentum can be expressed in terms of the 
eigenvalues of the Baxter operator as 
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Thus, the Schrodinger equation ( |1.2|) turns out to be equivalent to the eigenproblem for the 
Baxter Q— operator, Eq. ( p. 21) . 

The eigenvalues of the Baxter operator, Qq^qiu, u), have to fulfil the following three conditions: 

(i) Baxter equations: 

The function Qq^q{u, u) has to satisfy the holomorphic Baxter equation 



tN{u)Qq,q{u,u) = {u + is) Qq^q{u + i, u) + {u - is) Qq^q{u-i,u), (2.7) 

where 

t^iu) = 2u^ + q2U^-^ + ... + Qn (2.8) 

is the eigenvalue of the auxiliary transfer matrix with q = (g2, • • • , qN) denoting the eigenvalues 
of the holomorphic integrals of motion. The "lowest" integral of motion, q2, is related to the 
total SL{2, C) spin, h, of the system of particles 

q2 = -h{h-l) + Nsis-l), h = {l + nh)/2 + iUh (2.9) 

with Hh integer and Vh real. In addition, Qq^q{u, u) obeys the equation similar to ( p.7[ ) in the 
M— sector with s, h, and replaced, respectively, by 

-s = l-s\ h = l-h\ qk = ql (2.10) 

with k = 2, N. The function Qq q{u, u) does not depend on the total momentum of the state, 
P = iPiP) with p = —ij^k'^k P ~ "'^J^k^k^ invariant under h ^ 1 — h and 

h ^ 1 — h. Indeed, the Baxter operator Q{u,u) commutes with the generators of the SL{2,C) 
group and, as a consequence, its eigenvalue Qq^q{u, u) depends only on the SL(2, C) Casimir 
operators, like q2, which are symmetric under the above transformation of the total spin. 

{ii) Analytical properties: 

The spectral parameters satisfy ( ^.l] ) and their possible values can be parameterized as m = 
X — in/2 and u = X + in/2, with n arbitrary integer and A complex. Then, Qq^q{u, u) should be a 
meromorphic function of A with an infinite set of poles of the order not higher than situated 
at the points 

{m^ = ±i(s — m) , = ±i(s — rh)} , m,rh = 1,2, ... (2-11) 
The behaviour of Qq^q{u, u) in the vicinity of the pole at m = m = 1 can be parameterized as 



Qq,q{u^ + e, + e) = R^{q, q) 



1 , ^E^{q,q) 



(2.12) 



The functions R^{q,q) fix an overall normalization of the Baxter operator, while the residue 
functions E^[q,q) define the energy of the system (see Eqs. ( |2.15| ) and ( |2.18| ) below). 

{Hi) Asymptotic behaviour: 

In the above parameterization of the spectral parameters, Qq^q{u, u) should have the following 
asymptotic behaviour for | ImA| < 1/2 and ReA ^ oo 

Qq^q{\ - in/2, A + in/2) ~ e'®'^^^'^") xh+h-N{s+-s) ^ ^-^6,(5,5) ^i-h+i-h-Nis+s) ^ ^2.13) 
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with the phase 0/i(g, q) depending on the quantum numbers of the state and the total SL{2, C) 
spins h and h defined in and (|2.10|) . 



As we will show in Section 3, the Baxter equation ( p.7| ) supplemented with the additional 
conditions on the pole structure of its solutions, Eq. ( p.ll| ) , and asymptotic behaviour at infinity, 
Eq. ( |2.13| ), fixes uniquely the eigenvalues of the Baxter operator, Qq^q{u, u), and, therefore, allows 
us to determine the spectrum of the model. 

Additional properties of the function Qq^q{u, u) can be deduced from the symmetry of the 
model under permutations of particles. Apart from the cycle symmetry, Eq. ( p.5|) , the Hamilto- 
nian Ti^r is invariant under mirror permutations M \E'(ii, -Zat-i, z^q) = "^{zn, Z2, Zi) This 
transformation acts on the integrals of motion as — (— 1)*^^^ and maps the eigenstate of the 
Hamiltonian into another one with the same energy but different set of the quantum numbers 

[M^,,^] {zi,...,Zn-1,Zn) = ^-q,^gizi,...,ZN-l,ZN) (2.14) 

with — g = (g2, — 93, (— l)"'gn) and similar for q. For q2k+i = 0, or equivalently q = —q, 
Eq. ( p.l4| ) is replaced by M\l'g^g = (— l)^"'^~^"''\E'g^g. This property leads to the following parity 



relations for the residue functions R~^{q, q) defined in (p.l2| ) 

i?+(g, q)/R+{-q, -q) = e^'^^^t^^) (2.15) 

and for the eigenvalues of the Baxter operator 

QU-u, -u) = e'^^'i'^^ Q-,,-g{u, u) . (2.16) 

We recall that the spectral parameters u and u have to satisfy ( ^.Ij ). Examining the behaviour 
of the both sides of (|2.16|) around the pole at u = and u = uf and making use of Eq. ( |2.12|) 
one gets 

R^{q,q) = (-l)^e^^-(^'^")i?^(-g,-g), E^q,q) = -E^{-q,-q). (2.17) 

To obtain the expression for the energy E^i^q, q), we apply (^^) . Calculating the logarithmic 
derivative of Qq^q in the r.h.s. of ( p.4|) , we replace the function Qq,q{u ± i{l — s) , u ± i{l — s)) by 
its pole expansion (|2.12| ). Then, applying the second relation in (p.l7|) , one finds 



Er,{q, q) = E^{-q, -q) + (E+(g, q))* = Re [E+{-q, -q) + E+(g, q)] , (2.18) 

where the last relation follows from hermiticity of the Hamiltonian ( |1.5| ). We conclude from 
Eqs. (|2.18| ) and (|2.12| ), that in order to find the energy E^iq, q), one has to calculate the residue 



of Qq,q{u, u) at the (iV — l)th order pole at u = i{s — 1) and u = i{s — 1) . 

3. Quantization conditions 

Let us construct the solution to the Baxter equation ( |2.7| ) satisfying the additional conditions 



( p.l2|) and (|2.13|) . It proves convenient to use the following integral representation for Q 



q,q\"'7 



U,U] 



Qq,q{u,u) = [ ^^-»Z-»Q(Z,^), (3.1) 

J zz 

where integration goes over the two-dimensional i*— plane with z = z*. This ansatz is advanta- 
geous in many respects. Firstly, the condition (|2.1|) is automatically satisfied since it is only for 
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these values of the spectral parameters that the integral in the r.h.s. of (|3.1j ) is well-defined. 
Secondly, the functional Baxter equation on Qq^q{u, u) is translated into the A^th order differential 
equation for the function Q{z, z). Its derivation is based on the identity 

P{u) Qq,,iu, u) = ! — z-^^z-"^ Pi-izd,) Qiz, z) , (3.2) 
J zz 

with P{u) being a polynomial in u. Substituting ( |3.1| ) into the Baxter equation ( p.7| ) and applying 
one arrives at 



N 



z\zd,fz'-' + z-\zd^fz^-^ - 2{zd,f -Y,^'<ik{zd, 



k=2 



Q{z,z) = 0. (3.3) 



The ^—dependence of Q{z, z) is constrained by a similar equation in the antiholomorphic sector 
with s and qu replaced by s = 1 — s* and % = oX^ respectively. Finally, as we will show below, 
the remaining conditions on the analytical properties and asymptotic behaviour of Qq^q{u,u), 
Eqs. (|2.12| ) and (|2.13|) , become equivalent to a requirement for Q{z, z = z*) to be a single- valued 
function on the complex 2;— plane. 

Going through a standard analysis p4[, one finds that the differential equation ( |3.3| ) is of 



Fuchsian type with three regular singular points located at ^ = 0, ^ = 1 and z = 00. Defining 
linear independent solutions to Eq. (p.3|), Qa{z), and their antiholomorphic counterparts, Qh{z), 
we construct the general expression for the function Q{z, z) as 



TV 



Q{z,z)= Y,Qa{z)C,bQt{z), (3.4) 

a,b=l 

with Cab being an arbitrary mixing matrix. The functions Qa{z) and Qh{z) acquire a nontrivial 
monodromy around three singular points, z, z = 0, 1 and 00. For Q{z, z = z*) to be well-defined 
on the whole plane, the monodromy should cancel in the r.h.s. of ( p.4|) . This requirement leads 
to the set of nontrivial conditions on the matrix Cab- Solving them, we will be able not only to 
obtain the values of the mixing coefficients. Cab, but also determine the quantized values of the 
integrals of motion g^. 

We would like to point out that our approach to defining the function Q{z,z) in Eq. ( |3.4| ) is 
similar in many respects to the well-known approach to constructing the correlation functions in 
the minimal GET There, Q{z, z) plays the role of four-point correlation functions depending 
on the anharmonic ratios of the coordinates, z and z. The latter satisfy the differential equations 
("the null vector condition") similar to ( |3.3| ), in which the integrals of the motion are replaced 
by some combinations of the conformal weights of the primary fields. In the minimal GET, the 
conformal blocks Qa{z) and Qi,{z) are given by multiple contour integrals and their monodromy 
around the singular points z = 0, 1 and 00 can be found in a closed form. Going over to 
Eq. p.3| ), one finds (see Section 3.4 below) that similar representation exists only at = 2 
and it remains unclear whether it can be generalized for arbitrary number of particles N . Our 
subsequent analysis does not rely on such representation. 

To determine the function Q{z, z), we shall construct the r.h.s. of ( |3.4| ) in the vicinity of three 
singular points, z = 0, z = 1 and z = 00, and analytically continue the obtained expressions 
onto the whole 2;— plane. Additional simplification occurs due to the symmetry of the differential 
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equation ( p.3|) under the transformation 2; — > 1/z and Qk —>■ (— This symmetry is a mani- 
festation of a general property of the eigenvalues of the Baxter operator, Eq. (p.l6|) , which leads to 



Q,,,(^, z) = e^^-('''5) Q.,,.,{l/z, 1/z) . (3.5) 

Here, we indicated explicitly the dependence of the function Q{z, z) on the integrals of motion. 
Applying (|3.5| ), one can define z) around z = 00 from the solution at z = 0. 

3.1. Solution around z = 



Looking for the solution to ( p. 3D around z = in the form Q{z) ~ z*^, we find that the exponent 
a satisfies the indicial equation 

(a-l + s)^ = 0. (3.6) 

Since its solution, a = 1 — s, is A^— times degenerate, the small— 2; asymptotics of Q{z) contains 
terms ~ (In zY with k < N — 1. Let us define the fundamental set of linear independent solutions 
to (|3.3|) around 2; = as 



Qf\z) = z'-'m{z) 

m—l 

ur{z){\nzr-' + ^4_infc+i(z)(lnz)-'=-i 



k=l 



(3.7) 



with 2 < m < N and the binomial coefficients c'^_i = (m — l)!/(A;!(m — k — 1)!) inserted for later 
convenience. The functions Um{z) are given by power series 

00 

uUz) = l + J2z-u^r\<l), (3.8) 

n=l 

which converge uniformly inside the region \z\ < 1. Inserting ( p.7|) and ( p.8|) into ( p.3|) , one finds 
that the expansion coefficients (q) satisfy three-term (nonhomogeneous) recurrence relations 
with respect to n. To save space, we do not present here their explicit form. 

The fundamental set of solutions to the antiholomorphic differential equation, Q^^{z), can 
be obtained from ( p.7| ) by replacing s and by s = 1 — s* and (jk = respectively. Then, the 
general solution for Q{z, z) around 2 = is given by 

N 

Q{zrz)^'=' Yl Qt\z)&X4^{z). (3.9) 

m,rh=l 

The mixing matrix has to be chosen in such a way that Q{z^ z) should be single-valued 
at 2; = 0, or equivalently, the monodromy of Q^m {z) and Q^^{z = z*) around z = should 
cancel each other in the r.h.s. of ( p. 91) . According to ( p. 71) , the monodromy of Qm{z) is due to 
z^"*— factor and In 2;— terms. Taking into account that s — s = is an integer, we find that 



^The monodromy matrix, defined as Qn \z e'^'^'^) — MnkQl^\z), has a Jordan structure and, therefore, it can 
not be brought to a diagonal form upon redefinition of the fundamental basis. It is interesting to note that similar 
situation occurs in the operator algebra of primary fields in the so-called Logarithmic CFT [p6| . 
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the factor 



does not affect the monodromy of the r.h.s. of (|3.9| ). Then, for Q(z,z) to 



be single valued, it should depend only on \n(zz) rather than on \nz and \nz separately. It is 
straightforward to verify that this condition is satisfied provided that the matrix elements Cnm 
vanish below the main anti-diagonal, that is for n + m > + 1, and have the following form for 

n + m < N + 1 

N—n—m+l 



a 



(n- l)!(m- 1) 



E 

A;=0 



n+m—l 



(3.10) 



with (T, «!, otv-i being arbitrary complex parameters and oat = 1.0 

The mixing matrix depends on arbitrary complex parameters a and a^. We can 
fix the normalization of the function Q{z, z) by choosing the value of a. Since the resulting 
expression for the eigenvalue of the Baxter operator has to satisfy simultaneously two parity 
relations, Eqs. (|2.15|) and (|3.5| ), its value can not be arbitrary. As we will see in a moment, both 
relations are satisfied for a = exp^iO^^q, g)), with 6'jv(g, q) being the quasimomentum. Later, we 
will use (|0| ) to calculate the eigenvalues of ON{q,q) (see Eq. ( |3.24|) ). 

Substituting (|3.10|) into ( ^l9|) and taking into account ( ^^ , one gets small— z expansion of 
the function Q{z, z). The leading asymptotic behaviour for 2; — can be obtained by neglecting 
0{z) corrections in (|3.8| ). In this way, one finds 



\n^-'(zz) 



ln^-2(z^) 



_(A^-1)! (A^-2)! 



ln( zz) 

aN-i + ... H — ^2 + «! 



:i + 0{z,z)). 
(3.11) 

Using this expression, we can calculate the contribution of the small— ^ region to the eigenvalue 
of the Baxter operator, Eq. ( |3.1D . Introducing a cut-off, p <C 1, and replacing Q{z, z) in Eq. (|3.1|) 



by its expansion ( |3.11| ), one integrates term-by-term over the region \z\ < phj making use of the 
identity 



(Pz 



z-'^z-'''\n''(zz)z'' 



z\<p 



zz 



m—s—iu.m—s—iu 



-l)"n! 



lU 



\n+l 



+ C((m - s - m)°) 



with m and ffi positive integer. We find that, in agreement with the general properties of the 
Baxter operator, Eqs. ( p.ll| ) and ( p.l2|) , the function Qq^q{u,u) has poles of the order A^ at the 
points u = i{s — m) and u = i{s — fh). At m = m = 1 one finds from ( ^.11| ) 

g,,g-« + e, + e) = [I + le a^-i + ... + i^ef-^ ^2 + i^ef-' + 0(e^)] (3.12) 



with and defined in (|2.11|) . It is easy to see, using (|3.5|) and (|3.11|) , that the remaining 
poles of Qq^q{u,u) are located a.t u = —i{s — m) and u = —i{s — ffi) and they originate from 
integration in ( |3.1|) over the region of large \z\ > 1/p. 
Matching (^J^) into ( pJ2|) one gets 



R+{q,q) = - 



IT 

Jn 



-^i6»iv(<7,(j) 



E+(q,q) = aN-i{q, q) 



(3.13) 



We verify, using O^iq, q) = —6n{—(1i —q), that the obtained expression for R^{q, q) satisfies the 
parity relation (|2.15|) . According to their definition, Eq. (|3.10|) , the a— parameters are arbitrary 



•^To cancel the monodromy, it is enough to require that the sum over k should depend only on the sum n + m. 
We have chosen the sum in this particular form for later convenience (see Eq. (3.11)). 
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complex. We indicated in ( p.l3| ) the dependence of the a— parameters on the integrals of motion 
since we anticipate that their values will be fixed by the quantization conditions to be discussed 
below. Insertion of the second relation in ( |3.13| ) into ( |2.18| ) leads to the following remarkable 
expression for the energy 



EN{q, q) = Re [aN-i{-q, -q) + OLN-i{q, q)] 



(3.14) 



where q = {qk}, —q = {{~^)^qk} and similar for q. 

We conclude that the small— 2; asymptotics ( |3.11D leads to the correct analytical properties of 



the eigenvalues of the Baxter operator, Eq. (|2.11|) . Moreover, the energy of the system, E^iq, q) 



is related to the matrix elements of the mixing matrix ( 3.101) in the fundamental basis 



3.2. Solution around z = 1 

Substituting Q{z) ~ {z — 1)^ into ( p. 3D one obtains after some calculation the following indicial 
equation 

{h+l + h- Ns){h + 2-h- Ns)W(h-k)={), (3.15) 

k=0 

with the total SL{2, C) spin h defined in (|2.9| ). Since the solutions b = k with k = 0, N — 3 
differ from each other by an integer, one expects to encounter logarithmically enhanced terms 
~ ln(l — z). However, a close examination of ( |3.3|) reveals that the solutions to (|3.3|) do not 
contain such terms provided that h ^ {1 + nh)/2, or equivalently Im/i 7^ (see Eq. (p.9|)). At 
h = {1 + nh)/2 and Ims 7^ 0, the additional degeneracy occurs between the solutions to (|3.15|) , 
b = Ns — h ~ 1 and b = Ns + h — 2. It leads to the appearance of the terms ~ ln(l — z) in the 
asymptotics of Q{z) for z 1. Obviously, similar relations hold in the ^—sector. 

The fundamental set of solutions to Eq. ( |3.3D around z = 1 is defined similarly to ( |3.7] ). For 
Imh ^ one gets 

^f^\,^ 1^,. ^ AT ^ I, 

Vi{z) 

(3.16) 

with m = 3, N. The functions Vi{z) {i = 1, 2) and Vm{z) are given by the power series 

00 00 

v,{z) = i + Y,i^-zrv^\q), vUz) = i+ Yl i^-^rv^rTKq), (3.17) 

n=l n=N—m+l 

which converge uniformly inside the region |1 — 2;| < 1. We notice that Q''2 \z) can be obtained 
from the function Q^^\z) by replacing h ^ 1 — h. Substituting ( p.l6|) into ( p.3|) , one finds that 





= /-^ 


:i 


— 


Ns-h-l 




= z^-' 


:i 


— 


Ns+h~2 




= z^-' 


:i 


— 





the expansion coefficients f « ^ and f i"^^ satisfy the A^— term homogenous recurrence relations with 
respect to the index n.0 As was already mentioned, aXh = {l+nh)/2 the solutions Qi\{z) become 



z 



was included in the r.h.s. of ( 3.16| ) and ( ^.18 ) to simplify the form of the recurrence relations. 



Without this factor, the recursion will involve + 1 terms. 
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degenerate and one of them, Qi\z) for rih > 0, has to be redefined to include the additional 
ln(l — z)— term 



/i=(l+n,,)/2 



(3.18) 



Here the function V2{z) is the same as before, while Vi{z) = YlkLo'^^z'' coefficients Vk 

satisfy the A^— term recurrence relations with the boundary condition = 1. 



The fundamental set in the antiholomorphic sector, Q 



(1)/ 



is obtained from the functions 



Qn (z) by replacing s and h hy s = 1 — s* and h = 1 ~ h*, respectively. Among all functions 

: 1. As a 
1 can be 



in the fundamental set (|3.16|) only two, Q^i\z) and Q2'''{z), are not analytical at z 



consequence, a general solution for Q{z, z) possessing a trivial monodromy around z 
constructed as 



Q{z,z)'=' PnQ'^l\z)Q^^\z)+ [3,.uQ^i\z)Q^i\z)+ ^ g^^?(^)7™^Q 



7^(1) /^N 



N 



(3,19) 



m.m=3 



The /?— coefficients depend, in general, on the total spin h (and h = 1 — h*). They are chosen in 
( |3.19| ) in such a way that the symmetry of the eigenvalues of the Baxter operator under h ^ 1 — h 
is manifest. It is convenient to rewrite ( p. 191) in a matrix form as 



Ph 







Q(^,^)=^.C«.^, C«= I 

7 



(3.20) 



with 7 = 'jmm- The expansion (|3.19|) is valid only for Imh ^ 0. For h = {1 + nh)/2 the first two 
terms in the r.h.s. of ( p.l9| ) look differently in virtue of ( p.l8| ) 



Q{z,z) 



h={l+nh)/2 



A Q'i\^)Q^\z) + Q?i 



z)Q^^\-z) 



+ (32Q^\z)Q^^\z) + ... 



(3.21) 



where ellipses denote the remaining terms. Substituting (|3.19 ) into (3J.) and performing integra- 
tion over the region of |1 — 2;| <^ 1, one can find the asymptotic behaviour of Q(m, u) at large u 
and u. As we will show in Section 4, it turns out to be in agreement with the general properties 
of the Baxter operator, Eq. ( |2.13| ). 

The mixing matrix C'-^^ defined in (|3.20|) has a block- diagonal structure. It depends on 
2 + (A^ — 2)^ complex parameters [3h, Pi-h and which, in general, are some functions of the 
integrals of motion (g, q) to be fixed by the quantization conditions. Let us take into account that 
the function Q{z, z) has to satisfy the duality relation (|3.5| ). For |1 — -z] ^0, one can apply (|3.19|) 
to evaluate the both sides of ( |3.5| ) in terms of the mixing matrices C^^^(g, q) and C^^\—q, — g). 
This leads to the set of relations on the functions /3i(g, q) and jmrniQ, q)- 

To obtain these relations one uses the following property of the fundamental basis ( |3.16D 



N 



6=1 



Sab Q'b\z; q) 



(3.22) 



where Im(l/2;) > 0. Here, we indicated explicitly the dependence on the integrals of motion. 
Since the Q— functions in the both sides of ( |3.22| ) satisfy the same differential equation (p.3|), the 
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5*— matrix does not depend on z. As a consequence, one can evaluate the matrix elements Sah by 
examining the leading asymptotic behaviour of the both sides of ( p.22|) for 2; ^ 1. In this way, 
applying ( p.l6| ) and ( |3.17| ), one finds that the only nonvanishing matrix elements are given by 



-i-K^N s-h-\) 



22 



-iTT{Ns+h-2) 



k,k+m 



-1 



(A;-2s- 1) 



ml 



(3.23) 



with {x)m = r(x + m)/r(x), 3 < A; < and < m < N — k. Similar relations hold in the anti- 
holomorphic sector, Sn = e'^'^^'-'^-^l S22 = e'^(^^"+'^-2) and 5^^,^+^ = {-l)''-^{k - 2s - l)m/m\. 
Finally, we substitute (|3.19| ) and (|3.22|) into ( ^75|) and find 



P^{q,q) = e*«-('?''?)(-l)^-^+-'^/3,(-g,-g) 



N 



(3.24) 



n,n>3 



These relations imply that, similar to the energy, Eq. (|3.14|) , the eigenvalues of the quasimomen- 
tum, ON{q,q), can be obtained from the mixing matrix at 2; = 1. In particular, it follows from 
the first relation in ( p.24|) that the quasimomentum of the eigenstates with g2fc+i = ^2fc+i = 
{k = 1,2...), or equivalently q = —q, is equal to 



(3.25) 



since ph{q, q) = Ph{-q, -g). At = 2 one finds from (|]2|) that e^^^ = 



3.3. Transition matrices 



Eqs. ( ^.91 ) and ( |3.19| ) define the solution for Q{z, z) in the vicinity of 2; = and -2 = 1, respectively. 



To obtain the eigenvalues of the Baxter Q— operator, Eq. (|3.1|) , one has to sew (|3l9| ) and ( |3.19|) 
inside the region |1 — 2;| < 1, |2;| < 1 and, then, analytically continue the resulting expression 
for Q{z^ z) into the whole complex z— plane by making use of the duality relation ( p.5| ). As we 
will see in a moment, this can be done only for the special values of integrals of motion (g, q) 
satisfying the quantization conditions (see Eq. ( p.28| ) below). 

The sewing procedure is based on the relation between two fundamental sets of solutions, 
Eqs. ( |3.7| ) and ( ^3.161 ). Choosing z to be inside the region of convergence of the both series, 
Eqs. ( p.8|) and ( |3.17| ), we define the transition matrices fi(g) and Vt{q) 



N N 

Q^n\z) = Y.^nm{q)Q^\z) , 4n {z) = Y.^nm{q)4^ {z) • (3.26) 

m=l m=l 

Since the functions Qn\z) and Qm{z) satisfy the same differential equation ( |3.3|) , the transition 
matrices are independent. For the fundamental set of solutions, Eqs. ( |3.7| ) and ( ^.16| ), these 
matrices are uniquely fixed and they can be calculated as ||14| 

m = iyW[iy(i)]-i , = dlQ^\z,) (3.27) 

with i = 0,1 and Zq being some reference point, say Zq = 1/2, and similar for fi(g). The 
resulting expressions for the matrices ^l{q) and Q{q) take the form of infinite series in q and q, 
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respectively. The transition matrices in two sectors are related to each other as fl{q) = Qs{h,q) 
and = Qs{h,q). 

The transition matrices allow us to analytically continue the solutions ( |3.9| ) valid for |2;| < 1 
to the region |1 — 2:| < 1. Substituting (|3.26| ) into (|3.9|) and matching the result into (|3.19|) , 



we find that the two expressions for the function Q{z,z), Eqs. ( p^ ) and ( |3.19| ), can be sewed 



together provided that the mixing matrices C*-"^ and C*-^-* satisfy the following relation 

C^'\q,q) = [n{q)fC^^\q,q)n{q). (3.28) 

This matrix equation provides the quantization conditions for the integrals of motion of the 
model, qk and qk with k = 3, N. In addition, it allows us to determine the matrices C^^^ and 
C^^^ and, as a consequence, evaluate the eigenvalues of the Baxter Q— operator, Eq. ( p.l|) . Indeed, 
replacing in ( p. 281) the mixing matrices by their expressions, Eqs. (|3.9|) and (|3.19|) , we obtain the 



system of A^^ equations involving (A^ — 1) a— parameters inside the matrix C^^\ 2 + (A^ — 2)^ 
parameters /3i^2 and inside the matrix C^^\ as well as (A^ — 2) integrals of motion ^3, g^r 
(we recall that qk = ql)- Thus, the system ( p.28| ) is overdetermined. It allows us to determine all 
parameters including the quantized q and, in addition, it provides (2A^ — 3) nontrivial consistency 
conditions on the obtained solutions. Additional consistency conditions follow from ( p.24| ). 

The solutions to the quantization conditions ( |3.28|) for different number of particles A^ and 



the emerging properties of the spectrum of the model will be described in details in Section 5. 



3.4. Special case: N = 2 

As was already mentioned, the solution to the differential equation ( |3.3|) for A^ = 2 admit 
representation in the form of contour integrals and, as a result, the quantization conditions 
( p.28| ) can be solved exactly. At A^ = 2, after the change of variables z = {x — l)/x and 
Q{z) = [x{l - xy^-'" ^ 



^y{x), Eq. (p.3|) takes the form of the Legendre's differential equation ^4 

r ^ 1 

0. (3.29) 



d / \ di 1/1 \ 
— x(l - x)— + h{}i - 1) 



y[x) 



Its general solution is well-known as ?/(x) = (iw w''~^(w—l)''^^(w—x)~'^, where the integration 
contour Cy^ has to be chosen in such a way that the integrand resumes its original value after 
encircling Cw Then, two linear independent solutions to ( p.29| ) are given by the Legendre's 
functions of the first and second kind, P_/i(2x — 1) and Q_/i(2x — 1), respectively. Using the 
relation between these functions 



-7icot{nh) P-h{2x - 1) = Q-hi^x - 1) - Qh-i{2x - 1) 



(3.30) 



and going back to the z— representation, we choose the fundamental set of solutions to Eq. (|3.3| ) 
as Qs{z; h) and Qs{z; 1 — h), where the notation was introduced foif] 



Qsiz; h) 



l-s 



Q 



-h 



l + z 



(3.31) 
and 



The properties of the function Qs{z\ h), including its relation to the fundamental set 
( p.l6|) , can be found in Appendix ^ 

^Obviously, this definition is ambiguous. Instead of using the Q— functions, one may define the fundamental 
set entirely in terms of the P— functions (see Eq. (|A.9|) below). 
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Following (|3.4|) , we construct Q{z, z) as a bilinear combination of the functions Qs{z] h) and 
Qs[z; 1 — h) and their antiholomorphic counterparts, Qs{z; h) and Qs{z; 1 — h). Requiring Q{z, z) 
to have a trivial monodromy around z = 1 and taking into account that Qs{z; /i) ~ (1 — z)"^^'^'^ 
for z — > 1 (see Eq. 

(3.32) 



j2D), one gets 

Q{z, z) = Ch Qs{z; h)Qs{z] h) + Ci^h Qs{z] 1 - h)Qs{z; 1-h). 



To fix the coefficients Ch and ci-h, one examines the small— z asymptotics of (|3.32| ) with a help of 
Eq. ( |A.1| ) and requires that the terms ~ In 2; In ^ should cancel and the coefficients in front of In z 
and Inz should be the same. Applying the identity — h) — ip{h) = TTcot{7Th) one finds that 
the both conditions are fulfilled provided that Ch = —Ci^h and cot{nh) = cot(7r/;,). The second 
relation is automatically satisfied thanks to the property of the SL{2, C) spins of the principal 
series, h — h = rih with rih integer. Choosing 



Ch = — tan(7r/i)(— r 

TC 



rih 



we obtain from (|3.32| ) and ( |A.1| ) the small— z behaviour of the function Q{z, z) as 
Q{z, z) = z^-'z^-'i-l)""^ I ln{zz) + 2 Re [tP{h) + ^'(1 - h) - 2tP{l)] + 0{z, z) } 



(3.33) 



(3.34) 



Taking into account that the quasimomentum of the N = 2 states is equal to e*^^ = 

Eq. ( |3.25| ), we find that this relation is in agreement with ( |3.11| ). One determines the a— parameter 

by matching (|3|) into ( plj) 



ai{h) = 2Re[i){h) + ^{l - h) - 2i){l)] . 



Finally, applying ( p.l4|) we calculate the energy at = 2 as 

'1 + \nh\ 



E2{h, h) = 2ai{h) = 8Re 



The ground state corresponds to /i = /i = 1/2, or equivalently nh = Uh = 0, 

minE2(/i,/i) = -16 In 2. 



(3.35) 



(3.36) 



(3.37) 



and it defines the intercept of the BFKL Pomeron [0. 

The exact solution at = 2 is based on the properties of the Legendre functions, Eq. (|3.31 



Going over to the systems with the number of particles A^ > 3, one encounters the following 
difficulties. Firstly, representation for the solution to ( p.3[ ) in the form of contour integrals does 
not exist or, at least, it is not warranted. Secondly, for A^ > 3 the function Q{z,z) depends 
on the integrals of motion, q and q, whose values should be determined from the quantization 
conditions ( p.28| ). As we will show in Section 5, both problems can be solved by using the power 
series representation for the fundamental set of solutions, Eqs. ( p.7| ) and ( ^.16| ). 

In this Section, we have demonstrated that in order for the eigenvalues of the Baxter Q— operator 
to possess the prescribed properties, Eqs. ( |2.7| ) - ( p. 131 ), the integrals of motion, q, have to satisfy 
the quantization conditions (|3.28| ). In this case, one can construct the function Q{z,z) in the 
vicinity of the singular points, z = 0, z = 1 and z = 00, and analytically continue it onto the 
whole 2— plane with a help of the transition matrices ( |3.26 ). The spectrum of the model - the 
energy and the quasimomentum, can be obtained from the mixing matrices C^^^ and C^^\ which 
define the asymptotic behaviour of Q{z, z) around 2; = and z = 1, respectively. 
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4. Eigenvalues of the Baxter Q— operator 



In the previous Section, we established the quantization conditions for the integrals of motion 
q and obtained the expression for the energy Ejq. Let us now construct the corresponding 
eigenstates ^2, ^tv), Eqs. 



rj) and (U 



The analysis is based on the method of Separated Variables (SoV) developed by Sklyanin 
2^. It allows us to find the integral representation for the eigenstates of the model by going 
over to the representation of the separated coordinates x = (xi, ...,xn-i) [§| 



(4.1) 



with z = (zi, ijv). Here, U 



p,X^ 



transformation and ${g^q}(x 

expression for Uj^xi^) arbitrary was found in 
N = 2 and iV = 3 were obtained in 



is a kernel of the unitary operator corresponding to this 
is the wave function in the separated coordinates. The explicit 

. Independently, similar expressions at 



Remarkable property of the SoV representation is that $|g g}(5) is factorized into the product 
of the eigenvalues of the Baxter Q— operator depending on different separated coordinates 



N-l r 



(*{<?,5}(=2))*=e 



n 

k=l 



Vis + ixk)T{s - ixk) 



Til- s + ixk)T{l 



ixk) 



N 



(4.2) 



where the additional factor e'^^N{q,q)/'i jg needed to ensure ( p.l4| ). Contrary to the z = (2;, z)— co- 
ordinates, the possible values of the separated coordinates x^ = {xk,Xk) are quantized as fol- 
lows g, |T6| 



Xk = T^k ^ , = i/fc + — , (4.3) 

with Uk real and integer. Integration on the space of separated variables in Eq. (|4.1| ) implies 
summation over integer Uk and integration over continuous Uk 



N-l 



d 



N-l 



-> 1 2 



X 



Wk 



n E 

k=l \nj.=— CX3 



+ [nk 



Hi 



duk 



74. 



fi{x) 



27T- 



{N-l)\ 



i^n 



Xk 



Xj\ , 



(4.4) 



j,fc=i 

j>k 



where \xk 

Eqs. (|4.1j ) and ( [4.2|) allow us to calculate the eigenfunctions of the model in terms of the 
eigenvalues of the Baxter Q— operator. By the construction, the latter have poles specified in 
( p.ll|) . One verifies, however, that they lie outside the integration contour in ([4.4|) and the 
integral in ( ^.1| ) is well-defined. Still, one can make use of the pole structure of Qq,q{u, u) by 
closing the integration contour over Vk into the upper (or lower) half-plane and calculating the 
asymptotics of the wave function \l/p'^{g^q}(2) in the different regions of the space. 

According to ( p.l|) and ( |3.4| ), the eigenvalue of the Baxter Q— operator is given by the following 
two-dimensional integral 



d^z 

zz 



N 



(4.5) 



a, 6=1 



Similar integrals have already appeared in the calculation of correlation functions in two-dimen- 
sional CFT. Applying the results of [^], one can convert Qq^q{u,u) into a sum of products of 
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holomorphic and antiholomorphic contour integrals, j]jdz z~^'^'^Qa{z) and J^dz z~^~'^'^Qf,{z), 
respectively, with the integration contours C and C starting and ending at one of the singu- 
lar points z, z = 0, 1 and oo. These contour integrals define the set of 2A^— functions of the 
(anti) holomorphic spectral parameters u and u. In analogy with the CFT, we shall refer to 
them as the Q— blocks. The resulting expression for Qq^q{u,u), Eq. ( |4.5| ), is given by a bilinear 
combination of blocks belonging to two sectors. 

We would like to stress that, contrary to Qq^q{u,u), the definition of the Q— blocks is am- 
biguous. The eigenvalues of the Baxter operator, Qq^q{u,u), stay invariant if one replaces the 
Q— blocks by their linear combinations and redefines appropriately the expansion coefficients Cab 
in the r.h.s. of ( |4.5| ). Making use of this ambiguity, one may look for the definition of the blocks, 
for which the expression for Qq^q{u, u) looks particularly simple. In this Section, we shall present 
such a definition. We will demonstrate that the eigenvalues of the Baxter operator can be ex- 
pressed in terms of only two Q— blocks - one in each sector, defined below in Eq. ( [4.12|) . We will 
also show that the quantization conditions for the integrals of motion and the expressions for the 
energy and the quasimomentum, established in Section 3, can be reformulated in terms of the 
Q— blocks. 



4.1. Decomposition over the Baxter blocks 

To proceed with calculation of ( [4.5[ ), one has to specify the mixing matrix Cab, as well as the 
functions Qa{z) and Qf,{z). These functions have to be defined uniformly on the whole complex 

plane with the cuts and their bilinear combination, Eq. (|3.4| ), should match (|3l9|) and (|3.19|) 
for z — > and z ^ 1, respectively. 

Let us choose Cab to be the mixing matrix introduced in (|3.19|) and define the functions Qa{z) 
in the vicinity of 2; = 1 as 

Cab = Cj^Jiq, q) , Qa{z- q) Q«(^; q) . (4.6) 



Analytical continuation of Qaiz; q) to the region z and z ^ 00 can be obtained from (|3.26|) 
and (|]2g) as 



TV 



Qa{z;q) '=J2[n-\q)UQ!^\z;q), 

b=l 
N 

Qaiz; q) ™ [^(-^) S^ab Qf\^h; -q) , (4.7) 



6=1 



with the matrix 5* given by (|3.23| ). The functions Q}j{z; q) are defined similarly. 

The functions Qa{z) defined in this way possess a nontrivial monodromy aX z = 1. Encircling 
the point z = 1 on the complex 2;— plane in anticlockwise direction, one calculates from ( |3.16| ) 
the corresponding monodromy matrix as 

Qa{z) ^ MabQbiz) , M = diag(e2-(^^-^),e2-(^^+^), 1, l) . (4.8) 



Unity entries in this matrix correspond to (A^ — 2) functions in the fundamental set (|3.16|) 
analytical at z = 1. Following |28|, the two-dimensional integral in Eq. (|4.5|) can be evaluated as 
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(see Appendix |C| for details) 



Qq,q{u, U) 



if [(i-M-)c")]„.r 

a,b=l ^ 



dz 



^ dz ■ - 



(4.9) 



with M"^ = M according to ( |4.8| ). Here we tacitly assumed that the bilinear combination Eq. ( |3.4| ) 
is a single-valued function on the complex plane. This implies, in particular, that the integrals 
of motion q have to satisfy the quantization conditions (|3.28|) . 

As follows from their definition, Eq. ( |4.7| ), the functions Qa{z) satisfy the relation ( |3.22| ) 



N 



Qa{l/z; -q) = ^ Sab Qb{z; q) 



(4.10) 



6=1 



which holds for arbitrary z such that Im(l/z) > 0. Changing the integration variable in ( [4.9|) , 
^ — *• 1/z, and applying ( [4.10|) one gets 



QUu.u) = ^Y.[S\l-M)C^'^]^J ^z^-Q^{z;-q) [ Q^{z; q) . (4.11) 

We recall that the matrices S, M and C^^^ were defined before in Eqs. ( |3.23| ), ( [4. 8] ) and ( |3.19| ), 
respectively. Notice also that = S"^. Substituting the monodromy matrix ( [4.8|) into ([4.11|) we 
find that among A^^— terms in the r.h.s. of (|4.11| ) only two (with a = b = 1 and a = b = 2) provide 
a nonvanishing contribution to Qq^q{u,u). These two terms correspond to the Q— functions 



nonanalytical at 2; = 1, Qi{z) 



\Ns-h-l 



and Q2{z) 



l-z 



Let us introduce notation for the holomorphic and antiholomorphic Baxter blocks 



Q{u] h,q) 
Q{u; h,q) 



1 dz 



T{Ns - h) 
1 

V{Ns - h) Jo z 



z'"'Qi{z--q), 



'%-^^Q,[-z-q) 



(4.12) 



with the normalization factors chosen for later convenience (see Eq. ( [4.17] ) below). Taking into 
account Eqs. (|3.23|) and ( p.l9| ), one gets from ( [4.11|) 



Qq,q{u,u) 



TC 



T{Ns - h) 
_ V{l-Ns + h) 

V{Ns -l + h) 



+ 



r(2 ~Ns-h) 



Ph{q, q)Q{u]h,q)Q{u- h,q) 
f3i-hiq, q) Qiu; l-h,q) Q{u; l-h,q) 



(4.13) 



Here, we used the fact that the functions Qi{z; q) and Q2{,z] q) (as well as Qi{z] q) and Q2{,z] q)) 
can be obtained one from another through substitution h ^ 1 — h. 

We notice that in Eq. (|4.9|) the integration contours over z and z are different, whereas in 
the original two-dimensional integral ( [4. 5] ) the both variables appear on equal footing. Repeat- 
ing the calculation of ( [4.5|) and making use of the monodromy of the functions Qi,{z; q) in the 
antiholomorphic sector 



^ MabQbiz) 



M = diagf Q-^-^i.N-s--h)^^-2.i[N-s+-h)^ ^ 



(4.14) 
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one arrives at another (through equivalent) expression for the eigenvalues of the Baxter operator 



+ 



T{Ns - h) 
T{l-Ns + h) 

V{Ns -l + h) 



T{2-Ns-h) 



Ph{q, q) Q{-U] h, -q) Q{-U] h, -q) 
(3i-h{q, q) Q{-u; 1-h, -q) Q{-u; 1-h, -q) 



. (4.15) 



One can verify the equivalence of Eqs. ( [4.15|) and (|4.13| ), by substituting ([4.15|) into the l.h.s. of 
( p.l6| ) and taking into account the second relation in ( |3.24| ). 

Thus, for a given set of the integrals of motion, (g, q) , satisfying the quantization conditions 
( p.28| ), the eigenvalue of the Baxter Q— operator, Qq^q{u,u), is unique. It is expressed in terms 
of two chiral blocks introduced in (^4.12|) and is given by Eq. ( |4.13| ). 



4.2. Properties of the blocks 



Let us show that the blocks Q{u] h, q) and Q{u\ h, q), defined in ( |4.12|) , have the following prop- 
erties: 



[i) Q{u; h, q) satisfies the chiral Baxter equation ( |2.7| ). Q{u\ h, q) obeys similar equation in the 
antiholomorphic sector. 



[ii) Q{u;h,q) and Q{u]h,q) are meromorphic functions [^, ||, [T^ on the complex u— and 
M— planes, respectively, with the only poles (which are all of order not higher than N) 
located at the points m~ = —i{s — m) and = ^(s — m), with m and fh positive integer. 
The same property can be expressed in a concise form as 



Q(m; h, q) = r^(l - s + iu)f{u) , Q{u; h,q) = T^{l-s- iu)f{u) 
with f{u) and f{u) being entire functions. 



(4.16) 



{in) At large u and m, away from the poles ( |4.16| ), that is for Re(l— s+iw) > and Re(l— s— i-u) > 
0, the functions Q{u; h, q) and Q{u; h, q) have the asymptotic behaviour 



Q{u;h,q) 



[lU] 



-Ns+h 



[i + 0{i/u)] , Q{u-Kq) 



[l + 0{l/u)]. (4.17) 



By the definition, Eq. ( ^.12[ ), the blocks Q{u] h, q) and Q{u] h, q) are related to the same universal 
function calculated for different values of the parameters 



Q{u] h, q) = Qs{u] h, q) , Q{u] h, q) = Qs{-U] h, -q) 



(4.18) 



To verify these properties one uses the integral representation for the blocks, Eq. ( 4.12|) . 
Then, the first property follows from the fact that the functions Qi{z) and Qi{z) satisfy the 
differential equation (p. 31). As to the second property, the poles of Q{u; h, q) and Q{u] h, q) come 
from integration in (|4.12| ) over the region of small z and z. The leading asymptotic behaviour 
of the functions Qi{z) and Qi{z) in this region can be obtained from the first relation in (|4.7|) 
as Qa{z) ~ z^~^ ln^~^ z and Qb{z) ~ z}~^ ln^~^ z. Finally, to obtain the asymptotics at infinity, 
Eq. ( [4.17|) , one integrates in ( [4.12|) over the region z — > 1 and z —>■ 1 and makes use of Eqs. (|4.7|) 
and ( CT ) to replace Qi{z) ~ (1 - z)^''^'~^ and Q^{z) ~ (1 - z)^'-~^-\ 
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The above three properties uniquely specify the blocks for Imh ^ 0. As we have seen in 
Section 3.2 (see Eq. (|3.18|) ), at Im/i = one of the fundamental solutions to the differential 
equation ( pT3| ) has to be redefined in order to avoid a degeneracy. One encounters the same 
problem trying to define the block Q{u\ h, q) aXh = (l+n/i)/2. In this case, the linear combination 
Q{u] h, q) + cQ{u; 1 — h,q) satisfies the three conditions on the Q— block for arbitrary c and, as 
a consequence, the blocks Q{u; h, q) and Q{u] 1 — h,q) become degenerate. The expression for 
the blocks at /i = (1 + nh)/2 can be found in the Appendix 0. 

Using the definition of the blocks Q{u; h, q) and Q{u; h, q), one can establish different useful 
relations between them. In particular, as shown in the Appendix the blocks in each sector 
satisfy nontrivial Wronskian relations (Eqs. (p.4|) and (p.5| )). Moreover, the blocks in two sectors 
are related to each other as (see Eq. ([B.3|)) 



Q{u;h,q) 



r(l -s + iu) 
T(s + iu) 



N 



{Q{u*-l-h,q)y 



(4.19) 



At = 2 the eigenvalues of the Baxter Q— operator and the Q— blocks can be expressed 
in terms of the 3F2— hypergeometric series of a unit argument. As was shown in Section 3.4, at 
N = 2 the function Qi{z) entering (|4.12|) is equal (up to an overall normalization) to the Legendre 
function of the second kind. Substituting ( |3.31| ) into ( [4. 121) and using integral representation of 



the Legendre functions ||2^, one obtains after some algebra 
1 



Qs{u; h) 



1 — s + iu,l — s + iu,l — h 
1 + s — h + iu,2 — s + iu — h 



3-^2 



s + iu, 1 — s + iu, 1 — h 
1 + s — h + iu,2 — s — h + iu 



(4.20) 

where r[...] denotes the ratio of the products of the F— functions with the arguments listed in the 
upper and lower rows, respectively. Together with ( ^.13| ) and ( [A.7| ) this leads to the following 
expression for the eigenvalues of the Baxter operator at = 2 



Q{u,u) 



X 



= 2(-l 

r(2s- 



^'^ tan(7r/i) 
h) 



Lrfi 



l + h) 



2s - h) 



Q{u;l-h)Q{u]l-h) , (4.21 



where Q{u; h) = Qs{u; h) and Q{u; h) = Qs{—U', h) are (non-normalized) Q— blocks at A^ = 2.Q 

For A^ > 3 one can calculate the blocks by replacing the functions Qi{z; q) in (^4.12|) by their 
expressions in terms of the fundamental solutions around z = Q and z = 1, Eqs. ( |4.7| ) and (^]§), 
respectively. As shown in the Appendix P, this leads to two different series representations for 
the block Q{u; h, q), Eqs. ( |B.8D and ( [B.10| ), which are valid in the different regions on the complex 
plane. 

We are now in position to demonstrate that the eigenvalues of the Baxter operator constructed 
in this Section have correct asymptotic behaviour at infinity, Eq. ( |2.13| ). Substituting ( |4.17| ) into 
( [4.13|) , we verify that Qq^q{u,u) satisfies Eq. (|2.13|) with the phase Oh(q', q) given by 



.2j0h(g,5) 



j3h{q,q) T{2-Ns-h)T{Ns-h) 



(4.22) 



(3i_h{q. q) F(l -Ns + h)V{Ns -l + h)' 
In particular, for h = h = l/2+zz//j and Vh ^ ^ one has I3h{q, q) ^ l/vh leading to e^^^i/^C?.?) = _x. 



"Due to the additional fact or T'^{1 - h)/{2T{2 - 2h)) in the r.h.s. of (|A 
infinity that differs from (4.17) by the same factor. 



this expression has asymptotics at 
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4.3. Quantization conditions from the Q— blocks 



As we have seen in the previous Section, the eigenvalues of the Baxter Q— operator can be ex- 
pressed in terms of the Q— blocks satisfying the conditions In distinction with Qg^g{u, u), 
the Q— blocks can be constructed for arbitrary values of the integrals of motion q, q. In this Sec- 
tion, we will show that the requirement for Qq^q{u, u) to have correct analytical properties leads 
to the quantization conditions for the integrals of motion which are equivalent to ( |3.28|) . 



The general expression for the eigenvalue of the Baxter operator Q 



q,q 



in terms of the 



blocks Q{u] h, q) and Q{u] 1 — h,q) and their antiholomorphic counterparts looks like 



Q 



q,q 



Ch Q{u; h, q) Q{u; h, q) + Ci-h Q{u; 1 - h,q) Q{u] 1 - h,q) 



(4.23) 
u, u) defined 



with Ch (and Ci-h) being arbitrary function of h and the integrals of motion. Qq^q 
in this way is symmetric under h —>■ 1 — h and h —>■ 1 — h, it satisfies the Baxter equations 
in the holomorphic and antiholomorphic sectors, and its asymptotic behaviour at infinity is in 
agreement with Eqs. ( |2.13| ). Therefore, it remains to show that Qq,q{u, u) has the correct structure 
of the poles, Eq. (|2.11|) . To this end, one applies ( [4.19| ) (see also (p.3| )) and rewrites ([4.23 
two equivalent forms 



m 



Qq,q[U,U} 



lU] 



r(s — iu) 
r(l -s + iu) 
vis + iu) 



N 



N 



ChQiu; h,q){Q{u*; 1 - h,q)Y + Ci-hQ{u; 1 - h,q){Q{u*; 

ChQiu; h,q){Qiu*; 1 - + c,_hQiu; 1 - h,q){Qiu*; /i, g))*|(4.24) 



The analytical properties of Qq^q{u, u) are now manifest - the poles of Qq^q{u, u) in u and u are 
generated by the F— functions in the first and the second relation, respectively. One deduces 
from ( [4.24| ) that for arbitrary Ch and Ci_/i the function Qq^q{u, u) has the A^— th order poles at 
u = i{s — m) and, separately, at m = — i(s — m) with m, m = Z+. 



Let us now require that the analytical properties of ([4.24|) should match similar properties of 
the Baxter Q— operator, Eq. ( |2.11 ). We remind that the operator Q('u,'u) is well-defined only if 
the spectral parameters u and u satisfy (|2.1|) . Imposing this condition, we find that ( [4.24|) has 
the prescribed poles, ( |2.11| ), plus additional "spurious" A^— th order poles located at 



{u = z(s + m — 1) , u = i{s — m)} , 
{u = —i{s — m) , u = —i{s + rh — 1)} 



(4.25) 



with m,m = 1,2,.... Thus, the coefficients Ch (and ci-h) and the integrals of motions {q,q) 
have to be chosen in such a way that ( ^4.24|) has to have vanishing residues at the poles ([4.25|) . 
Introducing the functions 



$(e) 



Q{i{s + €);h,q) 



Me) = = 



Q{-i{s + e); h,q) 



Qit{s + e);l-h,q)' Q^^^^s + e);l -h,q) ^^'^^^ 

and examining the first and the second relation in ( [4.24|) for {u = i{s+'m — l + e),u = i{s — m + e)} 
and {u = —i{s—m+e), u = i{s+m—l+e)}, respectively, as e 0, one finds that this requirement 
leads to 



Ch $(m - 1 + e) + ci_ft (<I>(m - 1 + e))* 
Cft $(m - 1 + e) + ci-h {^{m - 1 + e))" 



O(e^) 



(4.27) 



21 



with m and rh being positive integer. It follows from the Wronskian relations (see Appendix B, 
Eq. (|B.6D ), that the infinite system of equations ( [4. 27]) becomes equivalent to a single condition 
at m = m = 1 



0{e 



Cl-h 



my 



Ch $(e) 



(4.28) 



with the functions ^(e) and <l>(e) defined in ( |4.26|) . 

One finds from ( |4.28| ) that Ch/ci^h is a pure phase. Its value can be obtained by matching 
(|]23|) into (|213| ) at large u and u with a help of ( ^T7| ), 



Cl-h 
Ch 



Then, recalling the definition of the function $(e), Eq. ([4.261), we obtain from (f4.28|) 



arg 



Q{i{s + e)-h,q) 
Q{i{s + e)-l-h,q)_ 



nh + 1 



and similar relation for the antiholomorphic block 



arg 



Q{-i{s + e)]h,q) 
\_Q{-i{s + e)-l-h,q)\ 



TC 



Uh + l 



+ 



Qh{q.q) + 0{e 



N\ 



(4.29) 



(4.30) 



(4.31) 



with I and £ integer, such that, in general, i ^ i. Applying ([4.19|) , one can express (|]31]) entirely 
in terms of the holomorphic blocks. The blocks entering the relations (|4.3CI| ) and ( [4.31| ) are given 
by (i3). 

Expanding the both sides of ( [4. 30] ) and ( [4. 311 ) in powers of e, one obtains the (overdetermined) 
system of 2A^ real quantization conditions on — 2 complex charges ga, ...,gAr and real phase 
0/i(g, g). We will verify in Section 5 that their solutions are consistent with the quantization 
conditions ( p.28|) . 



4.4. Energy spectrum from the Q— blocks 

Let us show that the energy, Ei^{q,q), and quasimomentum, 0]\[{q,q), admit a simple represen- 
tation in terms of the Q— blocks. To this end, we introduce new blocks 

Qo{u]q) = ahQ{u;h,q) + ai^hQ{u;l - h,q) , 

Qo{u]q) = ahQ{u;h,q) + ai^hQ{u;l-h,q) , (4.32) 

which are symmetric under h 1 — h and h ^ 1 — h. 

We require that Qo{u; q) and Qq{u; q) should have the same poles as the Q— blocks, Eq. ([4.16|) , 
but of the order not be higher than A^ — 1. Applying ( [4.19|) , it is straightforward to verify that 
the linear combination of the Q— blocks in the r.h.s. of ( [4.32| ) has a vanishing residues at the A^th 
pole at Ui = —i{s — 1) and = i{s — 1) provided that (up to an overall normalization) 

ah = tan{TTh)[Q{—is]h,q)y, = ta.n{7Th){Q{is] h,q))* . (4.33) 

Then, the residues at the remaining A^th order poles, {u^,u^}, vanish automatically, since 
otherwise Qo{u; q) and Qq{u] q) will not satisfy the Baxter equations (|2.7] ). 
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Let us now consider the following expressions 

> 1 N 



Qi{u]q) 
Qi{u;q) 



r 1 -s + iu) 



r(s + iu 
Til - s 



■-\-\ N 



T(s — iu) 



{Qoiu*; q))* = (ai„^)*Q(M; h, q) + {aj,yQ{u- 1 - h, q), 
{Qoiu*; q))* = {ai.tTQiu; K q) + (a,.)*Q(u; 1 - h, (4-34) 



where in the r.h.s. we used (|B.3|) . In distinction with (^4.19|) , this transformation maps Qo{u; q) 
into another block Qi{u] q) and similar for Qq{u] q). The ratio of the F— functions in the r.h.s. of 
( [4.34|) ensures that Qi{u] q) and Qi{u; q) have the same pole structure as the Q— block, Eq. (|4.16|) 
and, in addition, 

Qi(z(s + n) + e;g) ~ e, Qi(-i(s + n) + e; g) ~ e (4.35) 

for e — i> and n nonnegative integer. This property plays a crucial role in our analysis. 

Combining together ([4.32|) and ( [4. 341) , one can express the Q— blocks in terms of new blocks 
Qq{u; q) and Qi{u; q). Substitution of the resulting expressions into ( [4. 231) yields 



Qq,q{u, u) = Ag^g Qo{u] q)Qoiu; q) + Bg^gQi{u; q)Qi{u; q) . 



(4.36) 



Notice that the cross-terms, like Qo{u', q)Qi{u] q) , do not appear in this expression. They have 
nonvanishing residues at the iVth order spurious poles ( |4.25| ) and their appearance in the r.h.s. 
of (|4.36|) is protected by the quantization conditions ( |4.30|) and (|4.31|) . The coefficients Ag^g and 
Bq g can be expressed in terms of Ch, at and ajj, defined in Eqs. (|4.29| ) and ( |4.33|) , but we will 
not need these expressions for our purposes. 

To find the quasimomentum p.6| ), one has to evaluate ( [4.36| ) at u = ±is and u = ±is. Since 
the second term in ( |4.36|) vanishes due to ( [4.35| ), one gets 



fc'iv^g, g) = z In — : - + im : 



(4.37) 



The calculation of the energy (|2.3|) is more cumbersome. One finds from ([4.36|) and ( |4.35| ) that 
d 



du 



\\iQq^q{u -is,u- is) = (\nQo{-is; q))' + (ln(5o(-zs; g))' + AQ , 



u=0 



(4.38) 



where prime denotes the derivative with respect to the spectral parameter and 

^ Bg,gQi{-is)'^,i-ts) ^ 1 Q,{-ts)Qo{-t{s - 1 + e)) 
A,?Qo{-is)Qo{-is) '^0 ie Qo{-zs)Qi{-i{s - 1 + e)) 



(4.39) 



Here, the second relation follows from the requirement for ( [4.36|) to have vanishing residues at 
the "spurious" poles ( [4.25| ), or equivalently Qg,g{—i{s — 1 + e), —i{s + e)) to be finite as e — > 0. 
The residue of Qi{—i{s — 1 + e)) at the A^th order pole at e = can be obtained from the first 
relation in ( |4.34D . To calculate the residue of Qoi—i^s — 1 + e)) at the (A^ — l)th order pole at 
e = one applies the Wronskian relation 



Qq{u + i; q)Qi{u; q) - Qo{u; q)Qi{u + i;q) = const x 



r(m — s) 
T(iu + s) 



N 



(4.40) 
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which follows from (|B.4|) and the definition of the Qq— and Qi— blocks, Eqs. ([4.32|) and ([4.34|) , 
respectively. The normalization constant in the r.h.s. of ( [4. 401) can be obtained for u = —i{s + e) 
from the comparison of the residues of the both sides of ( [4.40| ) at the A^th order pole in e. In 
this way, one finds from ( [4.38 ) 



^lnQ„,( 



u 



IS, u 



IS 



iN 



?i=o 1 — 2s 

Repeating similar calculation of the logarithmic derivative of ( [4.36|) at u 



+ {\riQ,{-i-s-q))' -{\n[Q,{-i-s-q)r)' . (4.41) 

is and u = is we find 



^ In Qq^q{u + is,u + is) 



iN 



u=0 



1 - 2s 



^ + (\nQo{^s■q)y - (\n[Qo{^s■q)Y 



(4.42) 



and obtain the following expression for the 



Finally, we substitute the last two relations into 
energy 

EN{q,q) = -2lm{\nQo{is;q)y + 2lm{\nQo{-is;q)y + En , (4.43) 

where Sn = 2iVRe ['?/'(2s) + '?/'(2 — 2s) — 2'?/;(l)]. One can further simplify this expression by 
using the relation (\nQq^q{—is, —is))' = — (ln(5-q,-g(is, is))' that follows from (|2.16|) . Together 
with (14411 ) and (g^D, it leads to 



ENiq,q) = -2Im(lnQo(^s;g))'-2Im(lnQo(^s;-g))' + E 



(0) 
N ' 



(4.44) 



(0) 

N 



2A^ Re ['?/'(2s) + '?/'(l — 2s) — 2'?/'(l)]. The symmetry of the energy E]\f{q,q) 



where E[ 

Ei\i{—q, —q) becomes manifest in this form. 

The results obtained in Sections 3 and 4 provide the basis for calculating the energy spectrum 
of the model. We have demonstrated that 

• The eigenvalues of the Baxter operator, Qq^q{u,u), possessing the correct analytical prop- 
erties and asymptotic behaviour at infinity (see Section 2) can be constructed only for the 
special values of the integrals of motion q and q satisfying the quantization conditions, 
Eqs. (g), ( CT) and {J^^. 

• The functions Qq^qiu, u) can be decomposed over the chiral Q— blocks, which depend on the 
spectral parameters and the integrals of motion in the (anti)holomorphic sector, Eq. ( 4.13 ) 
and ( |4.36| ). In distinction with Qq^q{u, it), the definition of the Q— blocks is ambiguous and, 
therefore, they do not have any physical meaning per se. 

• Once the quantization conditions for the integrals of motion are fulfilled, the corresponding 
eigenvalues of the Baxter operator are uniquely fixed (up to an overall normalization). They 
allow us to calculate the energy spectrum of the model by using two different expressions 
Eqs. ( |3.14| ) and ( [4.44| ), which lead to the same value of En- 



We would like to mention that we disagree on these points with the approach of Ref. [IC 
which the Baxter equation for noncompact magnet of spin s = 1 has been investigated. 



m 



5. Energy spectrum 

In this Section we solve the quantization conditions (|3.28|) and describe the spectrum of the 
Schrodinger equation ( |1.2| ) for different number of particles N. 
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For arbitrary A^, the spectrum of the model - the energy, EN{q,q), and the corresponding 
eigenstates, \l/p>{g are uniquely specified by the total set of quantum numbers in the holo- 

morphic and antiholomorphic sectors, q = {q2,(l3, ■■■,<1n) and q = (^2, ^3, ^Af), respectively. 
Since qk = ql, the total number of independent complex valued quantum numbers is equal to 
(A^ — 1). One of them, q2, fixes the total SL{2, C) spin of the state {h, h) defined in (|2.9|). Ac- 
cording to Eq. (^]9D, the quantized values of /i, or equivalently g2, are parameterized by integer 
Uh and real Uh 

g2K, u,,) = -^-{^ + iv,}l + Ns{s - 1) . (5.1) 

At A^ = 2 this becomes the only quantum number parameterizing the spectrum of the model 
( p.36|) . For A^ > 3 one has to deal with a bigger set of the quantum numbers. 



To find the spectrum of the integrals of motion g3,...,gjv for A^ > 3, one has to solve the 
overdetermined system of the quantization conditions ( p.28|) . As was explained in Section 3, 



their solutions give us the expressions for the a—, /3— and 7— parameters entering the mixing 
matrices C^^^ and C*^^\ Eqs. (|3.9|) and ( |3.19|) . Then, the corresponding energy, EN{q,q), and the 



quasimomentum, ON{q,q), are calculated by inserting these expressions into ( |3.14[ ) and ( |3.24| ), 
respectively. In distinction with the N = 2 case, Eqs. ( p.36 ) and ( 3.35 ), the resulting expres- 



sions for the spectrum of the model can not be expressed for arbitrary A^ in a closed analytical 
form. Nevertheless, the results that we are going to present in this Section exhibit a remarkable 
regularity and they suggest that there should exist a WKB-like description of the spectrum. 
To analyze the quantization conditions ( p.28|) , one has to specify the values of the SL{2,C) 



spins (s,s) defined in (|1.4|) . Since our main motivation for studying the noncompact Heisen- 
berg spin magnets came from high-energy QCD, we shall fix them to be the same as for the 
A^— reggeized gluons compound state, Eqs. ( p..l|) and (|1.2| ), 

s = 0, s=l. (5.2) 

Another specifics of QCD is that, solving the Schrodinger equation ( |1.2| ), one is mainly interested 



in finding the ground state. It is this state that provides the dominant contribution to the 
asymptotic behaviour of the scattering amplitudes ( |1.1D at high-energy. Additional constraint 
on the solutions to ( |1.2|) is imposed in QCD by the Bose symmetry of the A^— reggeized gluon 
compound states. As we will argue below, this condition selects among all eigenstates only those 
with the quasimomentum exp^iO^iq, q)) = I. It is automatically fulfilled for the ground state for 
which one expects 9%°^^'^{q, q) = 0. 

To solve the quantization conditions ( p.28| ), one needs the expressions for the transition 
matrices fl{q) and fi(g), defined in ( |3.26|) . For N = 2 they can be determined exactly from 
the properties of the Legendre functions (see Appendix A). For A^ > 3 the calculation of these 
matrices is based on the relation (|3.27| ), which in turn relies on the power series representation 
of the fundamental solutions, Eqs. (|3.7|) and ( 3. 16 ). In general, the resulting expressions for the 



transition matrices take the form of an infinite series in q. Solving the quantization conditions, 
we truncated infinite series in Eqs. ( |3.8|) and ( p. 171) and retained a sufficiently large number of 



terms (rimax ~ 10^). This allowed us to calculate the energy spectrum numerically with a rather 
high accuracy (see Table |^ below). 

As was already mentioned, the system of A^^ quantization conditions ( |3.28| ) is overdetermined. 
Using the subset of (A^^ — 3A^ + 5) conditions, one can calculate the a—, (3— and 7— parameters 
in terms of the charges q and define the following test function 

/(g2,g3, ...,giv) = Tr(TTt) , T = C^'\q,q) - [mf C^'Kq, Q) m ■ (5-3) 
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Obviously, f{q) is a positive definite function on tlie (A^ — 1)— dimensional complex moduli space 
q = {q2,q3, ■■■,qN)- The solutions to the quantization conditions (|3.28|) correspond to points on 
this space, in which the test function vanishes, /(g) = 0. Since the dimension of the moduli 
space increases with A^, the problem of finding zeros of f{q) becomes very nontrivial at higher 
N. To solve it, we applied the algorithm developed in |^0[.Q It allowed us to identify the zeros of 
the test function ( ^.31) by the "steepest descent" method starting from some reference point on 
the moduli space. As yet another test of our approach, we verified that the obtained expressions 
for the integrals of motion satisfy the conditions ( |4.30|) and (|4.31| ). 



5.1. Fine structure of the spectrum 

Before summarizing the results of our calculations for A^ > 3, let us describe a general structure 
of the spectrum. We find that for given total SL{2, C) spin of the system, h = {1 + nh)/2 + iuh, 
the quantization conditions ( |3.28| ) provide us with an infinite number of discrete values of the 
integrals of motion. They can be parameterized by the set of integers as 

Qk = qk{^h;nh,{i}) , En = EN{i^h]nh,{i}) , {^} = (^i, ^2, ^2(7v-2)) (5.4) 

with k = 3,...,N. The explicit form of this dependence and interpretation of integers {£} will 
be given below. Eqs. (|5.1|) and (|5.4|) imply that, as a function of the total spin h, the quantized 
values of qk form the family of one- dimensional continuous nonintersecting trajectories in the 
(A^— 1)— dimensional space of g = (g2, qs, qN)- The "proper time" along each trajectory, Uh, is 
defined by the imaginary part of the total spin, Imh = u^, whereas the integers nh,ii, , ■■■■,(.2(n-2) 
specify different members of the family. Each trajectory in the g— space induces the corresponding 
trajectory for the energy = E^iq, g). We recall that, in contrast with the integrals of motion, 
g and g, the energy takes only real values. 

In what follows, we shall separate all eigenstates of the model into two groups according to 
the value of the "highest" charge q^'- qN ^ In = 0. The reason for this is that, as we will 
show in Section 5.5, there exists an intrinsic relation between the eigenstates with g^v = and 
the eigenstates of the system with the number of particles equal to A^ — 1 and qN-i 7^ 0. Namely, 
they share the same spectrum of the energy and the integrals of motion g^ (with k = 2, ...,N — 1) 
and the corresponding wave functions are related to each other in a simple way. 

For the moment, we shall exclude the eigenstates with gjv = from our consideration and 
return to them in Section 5.5. We find that among all eigenstates of the A^— particle system with 
gjv 7^ 0, the minimal energy occurs at i^h = = 0, or equivalently h = 1/2, 

""'^'^ = min E^(z/,; n^, {£}) = En{0; 0, {^ground}) (5.5) 

It belongs to the special trajectory {^ground}, to which we shall refer as the ground state trajectory. 
The energy along this trajectory is a continuous function of and it approaches its minimal 
value, at z//i = 0. In the vicinity of z//j = 0, one finds an accumulation of the energy levels 

ENii^h, 0, {^ground}) - i^^r'^"' = a^i^l + 0{vt) , (5.6) 

with (Tat being the diffusion coefficients. 

^This method was implemented in the form of the Fortran-90 code, which is available from the authors upon 
request. 
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The spectrum of quantized q2, ...,qN possesses the following symmetry 



Qk ^ i-irqk ^ ql (5.7) 

with k = 2, ...,N. Here, the first relation follows from invariance of the Hamiltonian under mirror 
permutations of particles, Eq. ( p.l4 ). The second relation is a consequence of the symmetry of 



the model at s = and s = 1 under interchange of the z— and sectors, or equivalently qk ^ qk 
and s(s — 1) s{s — 1). The relation (|5.7|) implies that if the quantization conditions (|3.28|) 



are satisfied at some point {q^} on the moduli space, then the same holds true at the points 
{{-!)%}, {ql} and 

As we will see in the next Section, the spectum of quantized charges q has a hidden structure, 
which can be revealed by examining the distribution of the quantized values of the "highest" 
charge q]^^ ■ Since q]^^ is a multi- valued function of complex g^r, each eigenstate of the model will 
be represented on the complex g^^— plane by different points. Together with the symmetry 
property of the spectrum, Eq. (IsTTI) , this leads to the following transformation on the moduli 
space 

q]l^ exp (nik/N) q]l^ , [Q<k<2N] (5.8) 

with k integer for odd and even for even A^. It maps one of the eigenstates into itself or into 
another one with the same energy. 

In the rest of this Section we shall present the results of our calculations of the spectrum of 
the model for the number of particles 3 < A^ < 8. 

5.2. Quantum numbers of the = 3 states 

At A^ = 3 the eigenstates depend on the quantum number ^3, which is an eigenvalue of the 
operator ^3 defined in ( p.. 6] ). Some of the eigenvalues of this operator have been already found in 



|T5|, |16] using different methods. The eigenstates found in [jl4| have pure imaginary values of 
quantized ^3 and their quasimomentum is equal to ^^3 = 0. We will demonstrate in this Section, 
that contrary to the statements made in the spectrum of the operator q^ is not exhausted 



by pure imaginary values, even in the sector with ^3 = 0. 

Solving the quantization conditions ( p.28| ) at A^ = 3, we reconstructed the full spectrum of 



quantized ^3. We found that apart from pure imaginary ^3 calculated in |T^, |T^, the spectrum 



also contains (an infinite number of) complex values of ^3, including pure real ones. Enumerating 
the quantized ^3 according to their absolute value starting from the smallest one, we notice that 

1/3 

1 53! grows cubically with its number. This suggests to describe the spectrum in terms of q^ 
rather than 03. To illustrate this point, we present in Figure |l| the results of our calculations of 
quantized ^3 for the total S'L(2,C) spin of the system /i = 1/2, or equivalently = = 0. 
Similar picture emerges at /i = (1 + nh)/2 + ivh for 'f^h = 1, 2, ... and Vh real. 

The spectrum of quantized ^3^'^, shown in Figure |l], is in agreement with the symmetry 
properties ( p.7| ) and (|578|). Defining the fundamental domain as < arg(g3^^) < 7r/3, we find 
that the whole spectrum of quantized ^3^^ can be obtained by applying the transformations ( |5.7|) 
and (p. 81) to the points belonging to this domain. 



It is difficult do not notice a remarkable regularity in the distribution of quantized charges 
in Figure 0. Apart from a few points close to the origin, the quantized values of ^3''^ are located 
at the vertices of the lattice built from equilateral triangles. As a consequence, they can be 
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Figure 1: The spectrum of quantized ql^^ for the system of = 3 particles. The total SL{2, C) 
spin of the system is equal to /i = 1/2. 



parameterized as 



1/3 



^N=3 



(5.9) 



where ii and £i are integers, such that their sum ii + ii is even. Here, A3 denotes the lat- 
tice spacing. Its value can be calculated from the leading-order WKB solution of the Baxter 



equation 31 



A, 





dx 




_4V3vr J 


-00 Vl - . 


271 



.395175... 



(5.10) 



1 /3 

Quantized occupy the whole complex plane except the interior of the disk of the radius A3 



^3^/^ I > A3 



(5.11) 



The comparison of ( |5.9| ) with the exact expressions for ^3 at /i = 1/2 is shown in Figure ^ and 
Table [|.| We find that the expression ( ^.9|) describes the excited eigenstates with a high accuracy. 
The agreement becomes less impressive for the eigenstates with smaller ^3. For instance, for the 
ground state with iq^ = 0.205258... and the first excited state with ^3 = 0.368293... the accuracy of 
( |5.9|) is ~ 16% and ~ 10%, respectively. Eq. (|5.9|) can be systematically improved by including 
subleading WKB corrections. Notice that the same WKB formula (|5.9|) is valid not only at 
h = 1/2 but also for arbitrary spin h. In the latter case, it describes correctly the excited states 
with 1^3^^ I ^ IqI^'^I- 

According to ( |5.9| ), the quantized values of ^3 are parameterized by a pair of integers ii and 
£2 which define the coordinates on the lattice shown in Figures |l] and 0. Eq. (|5.9| ) provides 

*To save space, we truncated in the Table the last few digits in the obtained numerical values of the charges 
and the energy. 
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Figure 2: Comparison of the exact spectrum of quantized q-^ at h = 1/2 (crosses) with the 
WKB expression (|5.9| ) (circles). 



(^1,4) 


(0,2) 


(2,2) 


(4,2) 


(6,2) 


(8,2) 


(10,2) 


exact ^1/3 


.590 i 


.358 + .621 z 


.749 + .649 z 


1.150 + .664i 


1.551 + .672i 


1.951 + . 676 z 




.684 i 


.395 + .684 z 


.790 + .684 i 


1.186 + .684i 


1.581 + .684i 


1.976 + .684 z 




-.2472 


-.6910 


-1.7080 


-2.5847 


-3.3073 


-3.9071 



Table 1: Comparison of the exact spectrum of ^3^'^ at h = 1/2 with the approximate WKB 



expression (|5.9|) . The last line defines the corresponding energy Es{i'h = 0; 



the WKB approximation to the function 53(1^/1; ^h, -^i, -^2) introduced in ( |5.4|) . Calculating the 
quasimomentum of the eigenstate specified by the pair of integers, ii and £2, one finds that it is 
given by 

27r 

03{h,i2) = Y^i (mod27r). (5.12) 

In particular, the quasimomentum vanishes for the eigenstates with ii = (mod 3). It is easy 
to see from ( ^.9|) that the corresponding take, in general, complex values. There are, however, 
special cases, like £1 = or £2 = 0, when ^3 becomes, respectively, pure imaginary or real. The 



former values have been previously found in [|I4|, p!5[] . 

We recall that the spectrum of ^3, shown in Figure ^ corresponds to the total spin h = 1/2, 
or equivalently Uh = i^h = in Eq. ( p.4|) . In general, quantized qs depend on integer Uh and 
continuous Uh- For simplicity, we present here our results only at n/i = 0, or equivalently h = 
1/2 + iuh- For Uh = 1,2,3,... the spectrum of ^3 exhibits a similar structure.^ 

We find that quantized ^3 = ^3(^/1, 0, £1, £2) are continuous functions of z/^. For different 
integers £1 and £2, the functions q3{i^h,0, £1, £2) define an infinite set of trajectories in the three- 

^For pure imaginary q^, the n^— dependence has been studied in [pO[. 
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dimensional Re(g3''^), Im(g3^'^)) moduli space. The trajectories cross the hyperplane Uh = 
at the points shown in Figure Q and go to infinity for ±00. To illustrate the properties 

of these trajectories, three representatives, corresponding to (£1,^2) = (0,2), (2, 2) and (4,2), 
are shown in Figure ^. The quasimomentum 6*3 takes a constant value along each trajectory, 
^3(^1, £2) = 0,47r/3 and 27r/3, respectively. 




Figure 3: The dependence of quantized qz{vh\ ^1, ^2) on the total spin h = 1/2 + iVh. Three curves 
correspond to the trajectories with (£1, £2) = (0, 2) , (2, 2) and (4, 2). 



Let us now consider the energy spectrum at = 3. Since the energy is a function of the 
total spin h and the charge ^3, each trajectory ^3 = qz{yh'AiA2) shown in Figure |^ induces the 
corresponding trajectory for the energy, £^3 = E^{i'h'AiA2) ■ Solving the quantization conditions 
( 3. 28 ) and applying ( 3.14 ) we obtain the energy spectrum shown in Figure 





-0.5 
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-^3/4 _1.5 
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Figure 4: The energy spectrum corresponding to three trajectories shown in Figure |^. The 
ground state is located on the (0, 2)— trajectory at I'h = nh = 0. 



We find that the energy is a continuous function of Uh along each (£1, £2)— trajectory and it 
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approaches its minimal value, niin,y,^ E^iyh'-, ^1,^2), at Vh = 0, or equivalently h = 1/2. Examining 
the value of E^^O; ii, £2) for different sets of integers (^1,^2), we find that E^^O; ii, £2) increases 
as one goes towards larger \q^ {0;ii,i2)\ (see Table |I]). This implies that the ground state 
corresponds to the point (s) on the plane of quantized ql^^ (see Figure 0) closest to the origin.[^ 
It is easy to see from Figure |I]that, in total, there are six such points, (£1, £2) = (0, ±2), (±3, ±1) 
and (=f3, ±1). According to (|5.12| ), their quasimomentum is equal to zero. Going over from q^^^ 
to g3, we find that these six points define two nontrivial eigenstates, which have opposite values 
of and the same energy 

^^ground ^ ±0.205258... , = .988678... (5.13) 

This implies that at = 3 and qs ^ the ground state is double degenerate. As we will 
demonstrate below, this property is rather general - the ground state is double degenerate for 
the systems with odd number of particles A^ and q^ 7^ 0, but it is unique for even A^. Eq. ( |5.13| ) 



is in agreement with the results of the previous calculations [p!4| , [T^ . To verify that for 7^ 



the ground state occurs at Uh = 0, one has to compare ( p. 13 ) with the minimal energy in the 



sectors with higher Uh- At ri/j = 1 , 2 , 3 our results can be found in the first three columns of the 
Table § (see below). They indicate that for qs the minimal energy grows as the Lorentz spin 
of the system, Uh, increases. 

At A^ = 3 the ground state is located on the (0, 2)— trajectory at Uh = Uh = as shown in 
Figure ^. We notice that close to Uh = there is an accumulation of the energy levels. At small 
Pfi the energy of excited states is described by a general expression (|5.6| ) with the dispersion 
parameter (T3 given below in the Table |^. 

5.3. Quantum numbers of the A^^ = 4 states 

At A^ = 4 the spectrum of the model depends on two complex quantum numbers q^ and q^. Similar 



to the A^ = 3 case, we shall determine their spectrum by solving the quantization conditions ( |3.28|) 
under the additional condition ^4 7^ 0. The only difference with the A^ = 3 case is that one has 
to increase the dimension of the mixing matrices, C^^-* and C^^\ and recalculate the transition 
matrices, fl and fl, as functions of q^ and ^4. As before, to understand the structure of the 
spectrum at = 1/2, it becomes convenient to deal with multi- valued complex variables g^^. 
Then, each eigenstate will be represented on the complex g^^— plane by four different points. 

At h = 1/2 the quantized values of {q3,q4,) are parameterized by four real numbers and, 
as a consequence, they do not admit a simple pictorial representation. The results of our 
calculations are presented in Figure ^ There, each point on the complex g]^^— plane has 
additional complex coordinate defined by the quantized values of the charge ^3. The latter 
are not displayed for simphcity. It is convenient to separate the eigenstates into two sets 
according to their quasimomentum 64 = (27rA;)/4, k = 0,2 and = 1,3 as shown in Fig- 
ure 1^. We notice that some points have very close values of q^. Nevertheless, the charge 
53 and the energy E4 corresponding to these points are different. Here is an example of a 
pair of such states at h = 1/2 and 64 = 0: (^3 = 0,^4 = -2.185790,^4 = 12.563898) and 
(g3 = 1.524585i,g4 = —2.195368, i?4 = 12.383790). Aside from this spurious degeneracy, the 

^°As we will show below, this property is rather general and it holds for the energy as a function on the 
{N — 2)— dimensional space of the integrals of motion (ga, gjv)- 



31 



1- 




+ 


+ 




+ 
+ 




+ 


+ 


1- 


++ + + 
+ + + 


+ + 

++ # 


+ + 
+ + 


++ 
+ 


0.5- 




+ 


++ 








++ 


+ 


0.5- 


+ ^ 
+ ++ 




+H 


+ 
+ 







+ 


+ 








+ 


+ 

-H- 


0- 

1— 1 
1 — 1 


+ ++ 




+H 


+ 


0.5- 




+ 

+ 


++ 

+ 




+ 
+ 




++ 

+ 


+ 
+ 


-0.5- 


+ ^ 

+ + + 


# ++ 




+ 

+ 


-1- 


















-1- 


++ + + 


+ + 


+ + 


++ 




-1 




-0.5 









0.5 


1 




-1 -0.5 





0.5 


1 



Figure 5: The spectrum of the integrals of motion at N = 4 and the total spin h = 1/2. The 
left and the right panels correspond to the eigenstates with the quasimomentum e^^^ = ±1 and 
±i, respectively. 



spectrum of quantized g4 has many features in common with the spectrum of ga at = 3 shown 
in Figure |I]. 

We notice that quantized g4 are located close to the vertices of a square-like lattice. To verify 
this property we selected among all eigenstates of the = 4 system only those with h = 1/2, 
gs = and nonzero values of (see Figure |^). These states have the quasimomentum 6*4 = 
and they play an important role in our discussion as the ground state of the system has the same 
quantum numbers. The WKB analysis of the Baxter equation at = 4 leads to the following 
expression for quantized q^ at gs = (see Eq. ( |5.19|) below) 



[g4™(^ 



1,^2 



1/4 



AAr=4 ■ — + I — 



.V2 V2 

where the integers ii and £2 are such that their sum ii + £2 is even and the lattice spacing is 



(5.14) 



"43/4 






-1 Vl~x^\ 



20F 



0.423606... 



(5.15) 



As before, the leading-order WKB formula (|5.14 ) is valid only for \ql^^\ ^ 1^2''^ I- The comparison 



, 1/2 1 



of ( |5.14| ) with the exact results for ^4 at /i = 1/2 is shown in Figure ^ and Table ^ (see footnote ||). 
We find that quantized occupy the whole complex plane except the interior of the disk of 
radius A4 

|gri>A4, (5.16) 

and the WKB formula (^.141) describes their spectrum with a good accuracy. 

We recall that the points shown in Figure |^ describe the A^ = 4 eigenstates with h = 1/2, 
gs = and the quasimomentum 64^ = 0. Similar to the A^ = 3 case, the spectrum of g4 is 
parameterized by the pair of integers (£1, £2), Eq- ( I5.14D . Still, as one can see from Figure |^, there 
exist the eigenstates which have the same quasimomentum, close value of g4 but gs 7^ 0. Similar 
phenomenon also occurs for other values of the quasimomentum. In order to distinguish these 
additional eigenstates, one has to introduce the second pair of integers {£3, £4). The dependence 
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Figure 6: The spectrum of quantized for h = 1/2 and = 0. The exact results are shown by 
crosses and the WKB predictions based on Eq. ( p.l4[ ) are denoted by circles. 



(^1,4) 


(2,0) 


(2,2) 


(3,1) 


(4,0) 


(3,3) 


(4,2) 


(g^exact)l/4 


.626 


.520 + .520 i 


.847 + .268 i 


1.158 


.860 + .860 i 


1.159 + .574? 


(^WKB)l/4 


.599 


.599 + .599 i 


.899 + .299 i 


1.198 


.899 + .899 i 


1.198 + .599i 


-E4/4 


0.6742 


-1.3783 


-1.7919 


-2.8356 


-3.1410 


-3.3487 



Table 2: Comparison of the exact spectrum of at ^3 = and h = 1/2 with the approximate 
WKB expression (|5.14| ). The last line defines the exact energy £4^. 

of qs and 94 on (£3,^4) leads to the "fine splitting" of quantized q^ in Figure ^. As we will argue 
below, the WKB approach allows us to describe a nontrivial spectrum of ^3 and ^4. 

Let us now examine the energy spectrum at = 4. Calculating the energy of the eigenstates 
shown in Figure we find that i?4 > for all points on the g^^— plane except of four points with 
the coordinates {ii, £2) = (±2, 0) and (0, ±2). Due to a residual symmetry q]^'^ — > exp(zA;7r/2) g^^, 
Eq. ( [5.8| ), they describe a single eigenstate, which can be identified as the ground state of the 
N = 4 system with h = 1/2 and 

gf°-°d^O, gr""'' = 0.153589..., Ef°"°'' = -2.696640... (5.17) 

It has the quasimomentum 6'4 = and, in contrast with the = 3 case, it is unique. 

Going over from z//i = to nonzero values of u^, one finds that, similar to the A^ = 3 case (see 
Figure the eigenvalues of the integrals of motion fiow along the trajectories on the moduli 
space {vhiq^iqi}- Namely, each point shown in Figure |^ creates its own trajectory labelled by 
the set of integers {£}. Among all trajectories the one with (^1,^2) = (2,0) plays the special 
role as it contains the ground state of the model. We find that ^3 = Img4 = for arbitrary 
Vh on the ground state trajectory, whereas Re(g4) and vary with as shown in Figure |^. 
Accumulation of the energy levels next to the ground state at z//^ = is described by Eq. (|5.6|) 
with the dispersion parameter (74 given below in the Table |^. 
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^-2 2 4 -4 -2 2 

Figure 7: The dependence of the energy, —£'4/4, and the quantum number, ^4/^2, with 
1/4 + ul, on the total spin h = 1/2 + iuh along the ground state trajectory at = 4. 



The ground state ( |5.17| ) has the Lorentz spin nh = 0. In analogy with the = 3 case, we 



calculated the minimal energy in the sectors with n/^ = 1 , 2 , 3 (see the first four columns in the 
Table || below) and verified that it increases with Uh- 

As we have seen before, the spectrum of quantized charges at A^ = 3 and A^ = 4 is described 
by the simple formulae, Eqs. ( |5.9| ) and ( |5.14| ). Their derivation is based on the WKB approach 
to the eigenproblem for the Baxter operator [^. In this approach, Qq^q{x,x) is constructed as 
a quasiclassical wave function in the separated coordinates x = {x,x). The underlying classical 
dynamics is described by the spectral curve ( "equal energy level" equation) 

r^: i/2 = 4(x)-4x2^ = (4x^ + g2a:^-' + ... + giv)(g2a;'^-' + ... + giv), (5.18) 



where t^ix) was defined in ( |2.8| ) and y{x) = 2x^ sinhp^: is related to the holomorphic part of the 
momentum of a particle in the separated coordinates, Px- For arbitrary complex x, the equation 
( [5 .181 ) has two solutions for y{x). As a consequence, y{x) becomes a single- valued function on the 



Riemann surface corresponding to the complex curve F^y- This surface is constructed by gluing 
together two copies of the complex x— plane along the cuts [(Ji,(T2], [(J2n-3,o'2n-2\ running 
between the branching points (jj of the curve ( ^.18| ). The latter are defined as simple roots of 
the equation t\i{aj) = 4:a'j^ . Their positions on the complex plane depend on the values of the 
integrals of motion q2, qs, qN- In general, the Riemann surface defined in this way has a genus 
g = N — 2, which depends on the number of reggeons, A^. It is a sphere at A^ = 2, a torus at 
A^ = 3 and so on. 

Let us define on F^v the set of oriented closed a— and /?— cycles. The cycles encircle the 
cuts [cr2fc_i, a2k] with k = 1, ...,N — 2 and belong to the both sheets of F^r. The cycles Pk run from 
the cut [o"2Ar-3, cr2N-2] to [(T2k-i, o"2fc] ou the upper sheet and, then, back on the lower sheet. Then, 
classical trajectories of a particle in the separated coordinates correspond to wrapping around 
a— and cycles on the Riemann surface (|5.18| ). Requiring Q'^^^^x, x = x*) to be single- valued 
on the complex x— plane, one obtains the following WKB quantization conditions^ 

Re i dSo = Til2k-i , Re I dSo = 7r^2fc , (A; = 1, A^ - 2) (5.19) 



"^We are most grateful to A. Gorsky for eoUaboration on this point. 
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where (ik are integer and (ISq is the "action" differential on the curve (|5.18| ) |pT | 



Solving ( [5.19D , one can find the explicit expressions for the integrals of motion, Eq. ( |5.4| ). At 
= 3 and = 4, for \q]^^ \ ^ qI^"^, one arrives at ( |5.9| ) and ( |5.14| ). The general analysis of 



( |5.19 ) turns out to be rather involved and it will be presented in the forthcoming publication. 



5.4. Quantum numbers of the states with higher A^^ 

In this Section we will describe the results of our calculations of the spectrum of the model for 
higher A^ and ^ 0. As we have seen in the previous Sections, the structure of the spectrum 
gets more complicated as one increases A^. Therefore, instead of presenting a detailed description 
of the whole spectrum, as was done before for A^ = 3 and A^ = 4, we will restrict our analysis to 
the properties of the ground state trajectory only. 

Solving the quantization conditions (|3.28|) for A^ > 5, we found that the ground state trajec- 
tories have different properties for even and odd number of particles. Namely, for even A^ the 
integrals of motion with odd indices k vanish and for even k they take pure real values 

93 = Q's = ••• = , Img4 = Imge = ••• = 0, [A^ = even] (5.21) 

For odd A^, the integrals of motion q^ take nonvanishing values. They are pure imaginary for 
odd k and real otherwise 

Regs = Regs = ••• = 0, Img4 = Imge = ••• = , [A^ = odd] (5.22) 

We recall that the total spin of the system on the ground state trajectory is equal to h = 1/2+iuh, 
so that the quantized charges and the energy Ej^ are functions of Vh- At A^ = 3 and A^ = 4 
the z//i— dependence of gTv and E]^ is shown in Figures ^, || and |^, respectively. 

It follows from Eqs. ( ^.21| ) and ( ^.22|) that for even A^ the ground state trajectory is invari- 



ant under the symmetry transformations ( |5.7| ), whereas for odd A^ it is mapped into another 
trajectory with the different quantum numbers (— l)'^gfc but the same energy E^. This implies, 
that the ground state of the system is double degenerate for odd number of particles A^ and it 
is unique for even A^. The degeneracy is related to the properties of the system under mirror 
permutations of particles |^. 

The results of our calculations of the ground state for the system with the number of particles 
N < 8 are summarized in the Table ^ and Figure ||. We recall that for odd A^ there exists the 
second eigenstate with the same energy En and the charges (— l)'^gA:. 

The energy of the ground state En, as a function of the number of particles A^, has a number 
of interesting properties. E^ changes a sign as one increases the number of particles - it is 
negative for even A^ and positive for odd A^. Our results also indicate that the absolute value of 
the energy decreases with A^, \En\ ~ 1/A^ for N ^ oo (see Figure For large even and odd 
number of particles A^ it approaches the same asymptotic value E200 = -E200+1 = 0. 

These properties are unique (and quite unexpected) features of noncompact SL{2, C) Heisen- 
berg spin magnets. It is instructive to compare them with similar properties of the 5'L(2,R) 



Heisenberg magnets studied earlier in |32, B3i p3, Rq] in the relation with the Evolution Equa- 



tions for high-twist operators in high-energy QCD. There, the number of sites of the spin magnet 
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94 










-^7v/4 




= 


2 














2.772589 


16.8288 


= 


3 


.205258 












-.247170 


.9082 


= 


4 





.153589 










.674160 


1.3176 


= 


5 


.267682 


.039452 


.060243 








-.127516 


.4928 


A^ = 


6 





.281825 





.070488 






.394582 


.5644 


N = 


7 


.313072 


.070993 


.128455 


.008494 


.019502 




-.081410 


.3194 


N = 


8 





.391171 





.179077 





.030428 


.280987 


.3409 



Table 3: The exact quantum numbers, g, and the energy, E'at, of the ground state of A^ reggeized 
gluons in multi-colour QCD. The dispersion parameter, ctn-, defines the energy of the lowest 
excited states, Eq. (p.6|). 



is equal to the twist of the operators and the ground state energy defines the minimal anomalous 
dimension of these operators. Roughly speaking, the Hamiltonian of the S'L(2,M) magnet is 
given by the holomorphic part of the reggeon Hamiltonian ( |1.5[ ) with the only difference that the 
holomorphic z— coordinates take real values and the eigenstates \E'(zi, z^) are polynomials in 
z. One finds |]33|, ^ |3^, that for the S'L(2,M) spin magnet all charges gs, ...,qN take real 
quantized values, the ground state energy Ej^ is positive and it monotonically increases with 
A^. Thus, the properties of the SL{2, C) and SL{2, M) Heisenberg spin magnets turn out to be 
completely different. 



-En/4: 
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Figure 8: The dependence of the ground state energy, —E^/A, on the number of particles A^. 
The exact values of the energy are denoted by crosses. The upper and the lower dashed curves 
stand for the functions 1.8402/(A^ - 1.3143) and -2.0594/(A^ - 1.0877), respectively. 



The spectrum of the SL{2, C) magnet is gapless for arbitrary A^. Accumulation of the energy 
levels next to the ground state is described by (|5.6|) with the dispersion parameter given in 
the Table 3. The states belonging to the ground state trajectory have the total SL{2,C) spin 
h = 1/2 + iuh and their quantum numbers ga, ...,gAr satisfy ( |5.21| ) and ( [5.22| ) for even and odd 
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N, respectively. The explicit expressions for can be obtained in the WKB approach from 
Eq. (IH). 




Figure 9: The dependence of the energy —Ej^iyh) 1'^ and the "highest" integral of motion \qN\/<i2 
with q2 = (1/4 + z/^) on the total spin h = 1/2 + ivh along the ground state trajectory for different 
number of particles 2 < < 8. At large z//j, —E^ > ... > —E^ > —E2 on the left panel and 
\Q8/Q2\ < ■■■ < \(l3/(l2\ on the right panel. 



Another unusual feature of the SL{2, C) magnet can be revealed by examining the dependence 
En = E]\r{h'h) on the total SL{2,C) spin h = 1/2 + iu^ along the ground state trajectory for 
different number of particles A^. For 2 < A^ < 8 this dependence is shown in Figure ^. We notice 
that the flow of the energy E^ with Uh is such that the hierarchy of the energy levels at z//j = 
and large Uh is completely different. At z//i = the values of En{0) coincide with those depicted 
by crosses in Figure |[ For large Uh the system approaches a quasiclassical regime 0, in 
which the energy EN{vh) and the quantum numbers qn have a universal scaling behaviour 

E^(z//,) ~41n|g^|, |g^|~C^i/^ (5.23) 

with C ~ 0.52. As can be seen from the right panel in Figures ^ and this regime starts already 
at Ph ~ 2. These results suggest that the ground state of the SL{2, C) magnet has the properties 
of quantum antiferromagnet, whereas the excited states are ferromagnetic. 

5.5. Descendent states 

Until now, we have excluded from consideration the eigenstates with = 0. Let us now 
demonstrate that these states can be expressed in terms of the (A^ — 1)— particle eigenstates. 
That is why we shall refer to them as descendent states. 

To see this, one examines the Baxter equation ( p.TD for the SL{2, C) spins s = and s = 1. 
One finds from (|2.8|) that for Qn = Qn = the transfer matrix factorizes as tN{u) = utN^i{u) 
and the holomorphic Baxter equation ( |2.7| ) at s = takes the form 

tN^i{u) Qn{u, u) = u^~^ [Qn{u + i,u) + Qn{u - i, u)] . (5.24) 

Here, we introduced a subscript to indicate that Qn{u, u) corresponds to the A^— particle system. 
Eq. ( |5.24| ) coincides with the Baxter equation for the eigenvalue of the Baxter operator for the 
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system of {N — 1)— particles, Qn-i{u,u). The same happens in the antiholomorphic sector. One 
verifies that for s = 1 the function uQn{u, u) satisfies the same antiholomorphic Baxter equation 
as Qj\f^i{u,u). This suggests that up to a unessential normalization factor 

Q^^''='\u,u)=Qn-i{u,u)/u. (5.25) 

However, in order to identify the r.h.s. of this relation as the eigenvalue of the Baxter operator 
for qN = Qn = 0, one has to verify that Qn{u,u) has appropriate pole structure, Eqs. ( |2.12| ) 



and (|2.11|) . By the construction, the function Q]\[_i{u,u) possesses correct analytical properties. 
Qn{u,u) inherits (A^ — l)th order poles of the function Qn-i{u,u), and, in addition, it acquires 
a spurious pole a.t u = u = 0. Therefore, Eq. (|5.25| ) defines the eigenvalue of the Baxter operator 
only if the residue at this spurious pole vanishes, or equivalently Qn-i{0, 0) = 0. As we will show 
in a moment, this condition is satisfied provided that Qn-i and the quasimomentum of the 
{N — 1)— particle state defined by the function Qiy_i{u, u) is equal to e*^^-i = (—1)^. 

To evaluate QAr_i(0,0), we examine of the both sides of (5.24) at m = u = e as e — > 0. 



Taking into account that Qiv_i(u, u) is finite at m = m = and it has {N — l)th order poles at 
(m = ±«,m = 0), Eq. ( CT ), 

QN-l{±^ + e,e) = [1 + 0{e)] , (5.26) 



we find from (|5.24|) that 



g7v-iQ^-i(0,0) = i?+„i(g,g) + i?]^_i(g,g). (5.27) 

The value of the residue functions R^_^{q, q) depends on the normalization of the Baxter oper- 
ator, while their ratio is fixed by the quasimomentum of the state. Indeed, it follows from (|2.15| ) 



and ( gT^ ) that R^_^{q,q) / R'^_^{q,q) = (_l)^-i e^^^^-iW'^). Combining this relation together 



with (15.271 ), we obtain 

g^-iQjv-i(0,0) = [l-(-l)^e-^^^-i]i?+_i(g,g). (5.2^ 



This relation is rather general and it holds for all eigenstates of the (A^ — 1)— particle system 
with the single-particle SL{2, C) spin s = and s = 1. In particular, replacing A^ — A^ + 1 and 
putting = 0, we deduce from (]5.28|) that the descendent states have the quasimomentum 



J9n 



-1 



gjv=o 



(5.29) 



Going back to ( |5.28| ), one finds that for qN^i the condition Qn_i{0, 0) = is satisfied provided 
that the quasimomentum of the "ancestor" (A^ — 1)— particle state equals e*^^-i = (—1)^. 

Thus, starting from arbitrary (A^ — 1)— particle eigenstate with the quasimomentum e*^^-i = 
(— 1)^ and applying (|5.25|) , one can obtain the A^— particle solution to the Baxter equation with 
qN = 0. Remarkably enough, the A^— particle state defined in this way has exactly the same 
energy as its ancestor (A^ — 1)— particle state 



E]\r{q2, qN-i, qN = 0) = E]\r^i{q2, qN-i) 



(5.30) 
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Indeed, as follows from (|2.18| ) and (|2.12|) , at s = and s = 1 the energy En is related to the 
behaviour of u'^QN^u + i, u) around m = 0. According to ( |5.25|) , at ^at = this function coincides 
with u^^^Q^^iiu + i,u). 

We conclude that the spectrum of the A^— particle system contains (an infinite) number of 
descendent ^at = states, which have the same energy as (A^ — 1)— particle states. An example 
of such states at A^ = 3 and A^ = 4 can be found in the Tables § and ||. To establish the 
correspondence between two different sets of the states, it is convenient to define a new quantum 
number, "^n = {—l)^~^^e~'^^'^ , so that •t}'^ = 1. It follows from our analysis that there is the one-to- 
one correspondence between the (A^—1)— particle eigenstates with ^^at-i = 1 and Qn-i 7^ and the 
A^— particle eigenstates with = 1 and = 0.0 This agrees with the results found previously 
in [^. One can argue that all A^— particle eigenf unctions with the highest charge qn = can 
be obtained in this way. We demonstrated that the (A^ — 1)— particle state with i^N-i = 1 and 
Qn-i 7^ can be transformed into the A^— particle state with = 0. Is it possible to continue 
this process and construct the (A^ + 1)— particle state with Qn-i and qn = Qn+i = 0? In this 
case, applying ( ^.25| ) one gets Qn+i{u,u) = Qn-i{u,u) /u^ . To compensate the spurious pole at 
u = u = Q one has to require that Qn+i{^, 0) = 0, or equivalently Q]si-i{e, e) ~ e^. Substituting 
this relation into the Baxter equation (|5.24|) ioi u = u = e and comparing the small e— behaviour 
of the both sides, one obtains the relation between the residue functions E^{q,q), Eq. ( |2.12| ), 



EN_i{q,q) = E^_^{q,q). Together with ( p.l7|) and (|2.18| ), this leads to E^-iiq. 



12, 



qN-i) 



0. 



Thus, a positive answer to the above question would imply that the energy of the (A^ — 1)— particle 
state should vanish for arbitrary values of z//j, that is along the whole trajectory on the moduli 
space. We do not find such states in the spectrum of the model. 



This implies that if one defines a linear operator A that maps the subspace Vj 



=1) 



of the 



eigenstates of the (A^ — 1)— particle system with qN-i and "^n-i 
states with qN = and {}n = 1, 



A : V, 



(l?iV_l=l) 



(^JV = 1) 



N-1 



1 into the A^— particle 



(5.31) 



then this operator is nilpotent |^. In addition, it has the following properties 

A2 = 0, a ■ 7^jv-i = Hiv • a , A-Pat^i = -Pjv- A = (-l)^A, (5.32) 
where the last relation follows from = PAr(— 1)^'*'"'^ with Pat being the generator of cyclic 



permutations of A^ particles, Eq. (2^). Then, one can show that the wave function of the 
A^— particle state with gAf = is related to its "ancestor" (A^ — 1)— particle state as 



(<?iv=0) 
N ' 



1 



A^n 



TV 



+ 



,Zn-1) 

1 

zn-i,n 



Hn-I "^N-lizi, . . . , ^Af-l) 



(5.33) 



with = Zk — Zn and Un = J2n=o Piv)" being the projector onto eigenstates with 

the quasimomentum (|5.29| ), PatHat = (— l)^^^PAf. The explicit expressions at A^ = 3 and A^ = 4 



^^We would like to stress that the charges qk (with k — 2, N) are continuous functions of the parameter Vh 
defining the total spin h = {1 + nh)/2 + ivh- The statement qn = means that qN{i^h) vanishes for all Vh, that 
is along the whole trajectory on the moduli space. If gjv vanishes only for some vu, the corresponding eigenstate 
is not descendant. Similarly, qn-i 7^ means that qN-i{vh) does not vanish except may be for some distinct Vh- 
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look as follows 



Z23Z31 Z31Z12 Z12Z23 

*f=°^ = ^^3{z^,Z2,Zs)-^M^2,Zs,Z,) 
Z34Z4I Z12Z41 

+ ^^*3(^3,^4,^l)-^^^3(^4,i'l,^2), (5.34) 
^12^23 ^23^34 

where the states \E'2 and \E'3 have the quasimomentum e*^^ = _x and e*^^ = 1, respectively. 
According to (p.25|) , e*^^ = (—1)"* and, therefore, the eigenstate v|/^'^='~°^ can be defined only for 



odd Lorentz spins n^. For > 4 the spin 71^ of the state '^^^'^ can be arbitrary integer. 
Let us transform the descendent states (|5.33|) into the momentum space using the relation 

ki...kN ■ ^N{ki, fc^v) = / n ^^^fc ^^^^1' ' (^-35) 

fc=l 

where k = [k^ + iky)/2 is holomorphic component of = {kx,ky) and the additional factor in 
the l.h.s. was introduced for later convenience. One finds from ( |5.34| ) after some algebra that up 
to a normalization factor 

3 1 1 2 2 3 

$r°^ = %|^$3(^4 + fc^l,4fc^3)-^^^^$3(4 + ^4,fc^l,4) 

4 1 3 4 

+ %±^$3(^2 + 44fci)-%i^$3(fci + fc2,44). (5.36) 



1.2 h.2 " L2 L2 

^3 rvi rv2 

;('?iv=0) 



Generalization to arbitrary is straightforward. The expression for coincides with the 

= solution found at A^ = 3 in Ref. and for arbitrary A^ in Ref. by different methods. 



One can verify Eq. ( [5.33| ) in a number of different ways. For instance, applying the operators 
( |1.6|) to the both sides of ( ^.33| ) one finds that the state \E'j^'^~°^ is annihilated by the operator 



qn and has the same spectrum (with k = 2, ...,N — 1) as the state "^n-i entering the r.h.s. 
of ( |5.33| ). It is more tedious to verify that ( |5.33| ) is in agreement with ( [5.25| ). To show this. 



one examines the action of the Baxter operator on the state ( ^.33|) , Q(n, u)\&^^ "'(^i, . . . , zn) 



replaces the Q— operator by its integral representation and uses the Feynman diagram technique 
to calculate the emerging two-dimensional integrals (see Ref. ^ for details). Finally, one can 
check by explicit calculation that P^v ^5^"^"°^ = (-1)^+1^5^^=°^ and A -^5^^=°^ = 0, in agreement 
with Eqs. ( |5.29| ) and ( p.32| ), respectively. 

Using the results of the previous Sections, one can apply Eqs. ( ^.30| ) and ( ^.33|) and reconstruct 
the spectrum of the descendent states for arbitrary A^. As before, we shall identify among these 
states the one with the minimal energy, E^^^'^\qN = 0). In virtue of (|5.3CI|) , this amounts to 



finding the minimal energy in the spectrum of the (A^ — 1)— particle system with Qn-i 7^ and 
the quasimomentum e*^^-i = (—1)^. Obviously, it can not be smaller then the energy of the 
ground state, Ef^'^^'^, whose quasimomentum is 6%'^i'^ = (see Table |). The results of our 
calculations at A^ = 3 and A^ = 4 are given in the Tables ^ and |^, respectively. As was already 
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E3 


i 


^3,6 


4 





.2052 i 


0.9884 









1 


.3315 


1.5368 


1 








2 


-.4017 + .42012 


4.6077 


1 






3 


1.1766 i 


6.9592 





8 






Table 4: The dependence of the minimal energy on the Lorentz spin nh aX N = 3. The total 
SL{2, C) spin equals h = {1 + nh)/2 for all states except Uh = 2 when h = 3/2 + .176 i. Integer i 
defines the quasimomentum, 6*3 = 27r£/3. The subscript (d) refers to the descendant states (see 



Section 5.5). The states with ga 7^ are degenerate with respect to (|5.7| ). 
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Qs 


94 


E, 
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93, d 


E4,d 
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.1535 


-2.6964 





.2052 i 


0.9884 
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1 











2 


1.3659 


8.1080 
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2 





-.2869 


2.6268 





1.0236 


6.8888 


2 


3 


-.6951 


-.6337 


5.8836 


1 


1.1766i 


6.9592 


2 



Table 5: The dependence of the minimal energy on the Lorentz spin Uh at N = 4. The total 
SL{2, C) spin equals h = (l + n/i)/2. The last three columns correspond to the descendant states. 
The quasimomentum is defined as 6*4 = (27r^)/4. The states are degenerate with respect to 



mentioned, at = 3 the descendant states can be constructed only for odd n/j. For Uh = 1 



and n/i = 3 they are descendants of the N = 2 states with the energy (|3.36|) . At = 4 the 
^4 = states are descendants of the A^ = 3 states with the quasimomentum e*^^ = 1. Notice 
that the A^ = 4 state with nh = 1 and gs = ^4 = is not descendant since 94 7^ for Uh ^ (see 
footnote |T^). 

In general, for even N one has e*^'^"^ = 1 and, therefore, the eigenstate with the mininal 
energy, \Efj^^~°\ is a descendant of the ground state, \E'fJ°™'^. As follows from our results (see 
Table ^ and Figure its energy is positive, i?^™(g7v = 0) = Ef^°^^'^ > 0. For odd N one has 
qiBn-i — _i Q^j^d^ therefore, this state can not be a descendant of the ground state. In this case, 
the minimal value of the energy is given by E'^^^lq^ = 0) = and the corresponding integrals of 
motion are equal to g2 = 93 = ••• = Qn-i = 0. This eigenstate is located on the trajectory on the 
moduli space at h = 1 + ivh and Vh = 0. Going over to z//^ 7^ we find that q2, 94, Qn-i 7^ 
while qs = q5 = ... = qN-2 = for arbitrary Uh- The energy E^^q^ = 0) is a continuous function 
of Uh along this trajectory and it takes its minimal value, E^^"^{qj^ = 0) = 0, at Uh = 0. We 
notice that it is smaller than the energy of the ground state on the subspace with q^ 7^ 0. 

At ^2 = ••• = 9Af = the eigenvalue of the Baxter Q— operator can be easily found from 
the integral representation ( |3.1| ). Solving the differential equations (|3.3|), one can construct a 
single- valued function Q{z, z) as 



Q(''=°)(z, z) = [co In^'-'izz) + ci In^^-^zz) + ... + c^-i] , (5.37) 



with Ci being arbitrary coefficients. To fix the coefficients, one substitutes ( ^.37|) into (|3.1|) and 



compares the asymptotics of the resulting expression for Q{u, u) at large u and u with Eq. ( |2.13|) 
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for /i = 1, s = and s = 1. One finds that ci 



Q 



(g=o), 



u 



u 



N 



u, u] 



Cn^i = and up to a normalization factor 

(5.38) 



with integer n defined in ( |2.lD . In this expression, the numerator compensates a spurious pole at 
u = u = Substituting (^.381) into (|2^) one calculates the corresponding energy, Ejq{q = 0) = 0. 



6. Summary 

In this paper, we have continued the study of noncompact Heisenberg SL{2, C) spin magnets 
initiated in |^. Having solved this model, we obtained for the first time a complete description 
of the spectrum of the multi-reggeon compound states in QCD at large N^. 

From point of view of integrable models, the results presented in this paper provide an 
exact solution of the spectral problem for completely integrable quantum mechanical model of 
interacting spinning particles in two-dimensional space. A unique feature of this model, leading to 
many unusual properties of the energy spectrum, is that its quantum space is infinite-dimensional 
for finite and conventional methods, like the Algebraic Bethe Ansatz, are not applicable. To 
overcome this problem, we applied the method of the Baxter Q— operator developed in application 
to the SL{2, C) spin magnets in Solving the Baxter equations, we were able to find the exact 
expressions for the eigenvalues of the Q— operator. They allowed us to establish the quantization 
conditions for the integrals of motion and, finally, reconstruct the spectrum of the model. 

From point of view of high-energy QCD, we calculated the spectrum of the colour- singlet 
compound states built from A^ reggeized gluons for A^ < 8 in the multi-colour limit, Nc —>■ oo. The 
obtained expressions allowed us to reveal the general properties of the spectrum for arbitrary A^. 
Our analysis was based on the identification of the A^— reggeized gluon states in multi-colour QCD 
as the ground states for the noncompact Heisenberg magnet of the length A^ and the SL{2, C) 
spins (s = 0,s = 1). The identification however is not straightforward. The contribution of the 
A^— reggeon states to the scattering amplitude ( p..l|) takes the form A = is'^j^^ias)^ An with 

^^(s,t) = I d2zoe^^-^^($a(^o)|e-^-<'°' 9^)"'|<I>5(0)) , (6.1) 

where dk = d/dzk, the rapidity Y = Ins plays the role of the "evolution time". The wave 
functions \^a{b){zo)) = ^a{b){^i—zo, zn—zq) describe the coupling of A^— gluons to the scattered 
particles. The "integration fixes the momentum transfer, t = —fp'. The operators l/<9^ stand 
for two-dimensional transverse propagators of A^— gluons and the scalar product is taken with 
respect to the SL{2, C) scalar product ( p..8|) . The Hamiltonian H^^^^^ = H^^^^^ + ... + H'j^^^^^ 
describes the interaction between A^ reggeized gluons in multi-colour QCD and it is given by the 
sum of the BFKL kernels [||, 0. Notice that it is different from the Hamiltonian of the magnet, 
Eq. (P). 

The Hamiltonians Ti.'^^^^ and H.^ act on different Hilbert spaces. For the magnet, it coincides 
with the representation space of the principal series of the SL{2, C) endowed with the scalar 
product (pr^). For the QCD Hamiltonian, the choice of the scalar product is dictated by physical 
requirements that the wave functions $a(b) have to be normalizable and the Hamiltonian 
has to be bounded from below. These conditions are satisfied if one normalizes the eigenstates 
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of7^r°' as 



N 

N 



|^(QCD)||2 
\^ N II 



J l[d'z,\d,...dr,^T°\z)\\ (6.2) 



k=l 

Under this choice of the normahzation condition, Tij^'^^' is related to the Hamiltonian of the 
SL{2, C) magnet of the spin (s = 0, s = 1) as 

^(QCD) ^ ^^Q^ ^^ Q^y^^^ (91. ..a^) , (6.3) 

where dk = d/dzk- One can verify using ( |1.5|) and (|1.3|) that this transformation changes the spin 
in the antiholomorphic sector from s = 1 to s = 0.0 Thus, on the space of functions normahzable 
with respect to (|6.2| ), two Hamiltonians have the same energy spectrum and their eigenstates are 
related as \E'7v,g(^) = di...d]\f'^'^^'^°\z). It remains unclear however whether physical solutions for 
7^^'^°'' can be constructed on a larger class of functions. 

The transformation (|6.3|) is not well defined on the subspace of zero modes of the operator 



di ...On- However, these modes do not contribute to (|6.1| ) due to gauge-invariance of the wave 
functions ^a{b){zi, z^)- This property can be expressed as |$a{6)) = di ... dNl'^aib)) with the 
states l^^aib)) normahzable with respect to the SL{2, C) scalar product ( |1.8| ). Then, substituting 
]3D into (|6.1|), one gets 



^;v(s,t) = f rf2zoe-°^^(vl/,(ib)|e-^=^«-/4|vl/,(0)) = 5^s-^=^-('''«")/^^ 

q,q 

(6.4) 

where the scalar product is taken with respect to ( |l.(j| ). In the second relation we decomposed the 
Hamiltonian over the complete set of its eigenstates ^I/p,{g,5} which have the total momentum p and 
the quantum numbers q = {q2, ...jQn)- The sum is dominated by the ground state contribution. 



Taking into account ( |5.6|) and performing integration over z//j one arrives at ( |1 . 1| ) . 

The wave functions "^a^izi, zn) are Bose symmetric - they are invariant under interchange 
of any pair of reggeon coordinates and their colour indices. Therefore the r.h.s. of (|6.4|) receives 
a nonzero contribution only from the states ^^p-^^g^g}, which have the same Bose properties. In the 
multi-colour limit, the Bose symmetry is reduced to the invariance of the wave function under 
the cyclic and mirror permutations of reggeons. Since the wave function of the A^— reggeon state 
is factorized, as Nc ^ oo, into a product of the colour tensor and the function of reggeon coordi- 
nates, the both factors have to possess the same parity under the cyclic and mirror permutations 
simultaneously, P = 1 and M = ±1. Since the operators of the corresponding transformations, 
P and M, do not commute with each other, this requirement leads to the selection rules on the 
eigenstates of the noncompact Heisenberg magnet. 

By the construction, the eigenstates \E'p>^{g^g}(zi, ...,zn) diagonalize the operator of cyclic per- 
mutations, Eq. (p.5|) . Then, making use of ( p. 141) , the eigenstates of the operator of mirror 
permutations, M\I'l^^ = ± \E'i^^-|, can be defined as 

= = I [%<?,5}(^) ± , (6-5) 

where q = {g^} and —q = {(— l)^gfc}. Although these states do not diagonalize the integrals 
of motion, {q,q}, they are the eigenstates of the Hamiltonian having the same energy Ei<^ = 



13 



This corresponds to switching from strength tensors to gauge potentials in the description of gluon correlations. 
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En{(1, q) = E]\[{—q, —q). Using ON^—q, —q) = —ON^q, q), we find from ( pTsD that the operator of 
cychc permutations acts on them as 

= q)) ■ 4tL-}(^) + ^sin(^^(g, q)) ■ ^^l.yiS) . (6.6) 

Thus, the eigenstates diagonahze the operators P and M simultaneously only at 

sin(6'Ar(g, g)) = 0. Together with (|2.5|) , this condition selects among all eigenstates of the mag- 
net only those with the quasimomentum e^^^^'^'^^ = ±1 for even and e^^^^'^i'^) = 1 for odd A^. 
Obviously, it is satisfied for the ground state, 6%"^^'^{q, q) = 0. 

We conclude that the ground state of the magnet has a definite parity with respect to the 
cyclic and mirror permutations simultaneously and, therefore, its wave function \E'^^-|(2) can 
be identified as the i*— dependent part of the full wave function of the compound states of N 
reggeized gluons a.s Nc ^ oo. The ground state of the magnet has different properties for even 
and odd number of particles. For even N, its quantum numbers satisfy q2k+i = 0, Eq (|5.21|) , 
and, as a consequence, two sets of the quantum numbers, (g, q) and (— g, —q) coincide leading 
to \E't^^-|(z) = ^^^q^q}{z) and \E'l|^-j(5) = 0. For odd A^, the ground state in the sector with 
gAT 7^ is double degenerate. The degeneracy occurs due to the symmetry of the energy E]\f 
under q2k+i — g2fc+i, Eq. ( p.l8|) . This allows us to construct two mutually orthogonal ground 



states, \E'L|^-|(z), which are invariant under the cyclic permutations, 9n = 0, and possess a 
definite parity under the mirror permutations, M \E'^^ -| (2) = ±\l/^^ In virtue of the Bose 
symmetry, the colour part of the wave function of the reggeized gluon state should have the 
same parity under the charge conjugation, t° — * —(t'^)'^, with t°- being the SU{Nc) generators 
in the quark representation. This allows us to distinguish the ground states according to their 
C— parity. For even A^ the ground states with the parity M = 1 have the same C— parity as 
the Pomeron, C = 1. For odd A^, in the sector with gA? 7^ 0, the ground states with the parity 
M = 1 and M = — 1 have the C— parity of the Odderon [1^, C = —1, and the Pomeron, C = 1, 
respectively. For the descendant states, g^r = 0, one deduces from ( ^.291) that the physical states 



with 6]s[ = can be constructed only for odd A^ [^. Their minimal energy is -E'^^in'' ~ ^ 



the corresponding state with g2 = ... = gA^ = has the quantum numbers of the Odderon. 

Our results indicate that, in the multi-colour limit, in the Pomeron sector, only compound 
states built from even number of reggeized gluons A^ provide the contribution to the scattering 
amplitude .4(s, t)/s, Eq. rising with the energy s. Their intercept = 1 — agE^/i is 

bigger than one, but it decreases at large A^asaAr — 1~1/A^. 

In the Odderon sector, the situation is different. Depending on the value of the "highest" 
charge, gAr 7^ and gAr = 0, one can construct two different solutions for the Odderon state. 
For the first solution, gAr 7^ 0, the intercept of the compound states is smaller than one, but 
it increases with N as — ^ ^ —'^/N. As a consequence, the contribution to the scattering 
amplitude ( |1.1| ) from the A^ = 3 state ("bare Odderon") is subdominant at high-energy with 
respect to the contribution of the A^ = 5 state and so on. The high-energy asymptotics of the 
scattering amplitude ( |1.1| ) in the Odderon sector with gA? 7^ is governed, as A'c — > 00, by the 
contribution of the states with an arbitrary large odd number of reggeized gluons. It increases 
the effective value of the Odderon intercept and leads to aoddcron = ^200+1 = 1- For the second 
solution, gAT = 0, the intercept of the compound states equals 1 for arbitrary odd A^. At A^ = 3 
such state was first constructed in To calculate the contribution of the gA? = states to the 



scattering amplitude, one has to resum in (1.1) an infinite number of terms with A^ = 3, 5 



They have the same energy behaviour ~ s^/(crAr In s)^/^, with the dispersion parameter ajy, which 
scales at large N as ~ 1/A^^ and, therefore, enhances the contribution of higher reggeon states. 
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Thus, the two solutions, g^v 7^ and qn = 0, lead to the same value of the Odderon intercept, 
ttodderon = 1, but the properties of the underlying Odderon states are quite different. The 



Odderon state with Qn ^ does not couple to a point-like hadronic impact factors of the form 

~ 5{zi — Z2)(f){z2, + [cyclic permutations], like the one for the 7* r/c transition, 

whereas the Odderon state with ^at = provides a nontrivial contribution |3^. Another difference 
comes from the analysis of the dependence of the scattering amplitude on the invariant mass of 
one of the scattered hadrons, Q^. One can show that in the limit xbj = Q^/s ^ 1 and — ^ 00, 
the Odderon states with giv 7^ and = provide a contribution to the scattering amplitude 
~ l/Q^ of the twist-4 {p = 4) and twist-3 (p = 3), respectively. Both properties have to do 
with the fact that the wave function of the A^— reggeon state, \E'p'^{q^g}(ii, ^at), vanishes as 
\zk — Zk+i\ — ^> for 7^ and it stays finite for = 0. It remains unclear which of these 
solutions corresponds to a physical Odderon state in QCD. 

We found that in the Pomeron and the Odderon sectors, the intercept of the A^— reggeon states 
approaches the same value aoo = 1 as N —>■ 00 and their contribution to the scattering amplitude, 
^(s, t)/s, ceases to depend on the energy s as A^ — *• cxo. It is interesting to notice that this result 
has been anticipated a long time ago within the bootstrap approach B. It is also in agreement 



with the upper bound on the energy of the compound reggeized gluon states established in |T2 



We would like to remind, however, that the calculations were performed in the multi-colour limit 
and the important question remains: may the nonplanar corrections change the A^— dependence 
of the energy of the compound reggeized gluon states? One expects that the nonplanar 
corrections to the reggeon Hamiltonian will break integrability of the Schrodinger equation (|1.2|) 
and calculations will be more involved. This problem deserves additional studies. 

Analyzing the high-energy asymptotics of the scattering amplitudes, one is trying to identify 
the effective theory, which describes the QCD dynamics in the Regge limit. The main objects of 
this effective theory are the A^ = 2, 3, ... reggeon compound states constructed in this paper. In 
the generalized leading logarithmic approximation, these states propagate between the scattered 
hadrons and do not interact with each other. In the topological l/A"^— expansion 0, these states 
emerge from the summation of cylinder-like diagrams, whose walls are built from the reggeized 
gluons. These diagrams can be interpreted as describing the propagation a closed string between 
two scattered hadrons. This suggests that the effective dynamics of the multi-reggeon compound 
states in multi-colour QCD has to admit a stringy representation EO]. 
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Note added 

After this paper has been submitted for publication, the preprint by H. de Vega and L. Lipatov 
(dVL) appeared, [arXiv: [hep-ph /0204245|| , in which the spectrum of the A^ = 3 and A^ = 4 
reggeon compound states was investigated. It contains a number of statements, which are in 
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contradiction with our results. We would like to comment on them below. 

In the dVL paper, it is claimed that, contrary to our findings (see Figures |l|, ^and Tables |l], ||), 
the quantized valued of the integrals of motion for arbitrary can take only pure imaginary 
values for odd k and real values for even k. In particular, at = 3 this would imply that 
Regs = for all eigenstates of the model. Moreover, for complex values of the charges found 
in our paper, the dVL approach leads to complex values of the energy Ej^. The authors attributed 
a disagreement between their and our results to the fact that the quantization procedure proposed 
in our paper is erroneous. They did not offer, however, any further explanations and refer instead 



to [H. In order to check the above assertions and, at the same time, to test the dVL approach. 



we decided to verify our results at = 3 using the approach by R. Janik and J. Wosiek |jT^. We 



found that all eigenstates at = 3 constructed in our paper, including those with Regs ^ 0, 
satisfy the quantization conditions formulated in ||14[-0 example, we present two such 



states with £2) = (3, 3) and (3, 5), which have the total spin h = 1/2 and the quasimomentum 
^3 = (see Eqs. (|;9|) and (CT)) 

^3(3, 3) = -1.475327 , = 8.469248 , 

^3(3, 5) = -4.752678 - 3.048722 i , = 13.850368 . 

In the notations of the wave functions of these states are specified by the parameters, 
corresp ondingly, 

0^^ = 0.4332, = -0.6345 - 0.0361 i, 7^^ = 0.6391, 

a,^ = 0.0895 + 0.0007i, = 0.8582 , 7^^ = 0.5055 + 0.0037i . 

Moreover, one can argue that the set of eigenstates with Regs = can not be complete. If it 
were complete, the operator gs + gs defined in Eqs. ( |1.6|) and (|1.7| ) would be identically equal 



to zero on the SL{2,C) representation space ( |1.8| ). Applying this operator to an arbitrary test 
function on this space "^{zi, Z2, z^), one verifies that (gs + gs)^ 7^ 0. In a similar manner, one 
can show that for higher A^ the eigenvalues of the integrals of motion gs, gTv can not take only 
real or pure imaginary values. 

Another issue concerns the ground state at A^ = 4, Eq. (|5.17|) . In the first version of the dVL 
paper, it was claimed that it is located aX h = 1 and its energy is smaller than the energy of the 
state defined in Eq. (|5.17| ) above. Later, in the second version of the paper, the authors found 
yet another state at h = 3/2, which has even smaller energy. We remind that for the principal 
series of the SL{2,C) group the total spin h has the form (^l9| ), so that h = 1/2,1 and 3/2 
correspond to z//i = and the Lorentz spin Uh = ,1 and 2, respectively. Thus, the results of the 
dVL paper imply that the minimal energy of the system of A^ = 4 particles decreases as the total 
angular momentum of their rotation on two-dimensional plane, Uh, increases. Our results (see 
Table ^, 4th column) indicate that the dependence is just opposite. In addition, in our approach 
we do not find such states with h = 1 and h = 3/2 in the spectrum of the SL{2, C) magnet and 
the ground state occurs at h = 1/2, Eq. ( p.l7| ). Analyzing the quantization conditions ( p.28| ), 
we were able to identify the physical meaning of the solutions found in the dVL paper. They 
fulfil the quantization conditions at A^ = 4, but have unusual quantum numbers, h = h = 1 and 
h = h = 3/2, which do not match the principal series of the SL{2, C) group, Eq. (|2.9| ). One can 
show that these states are located on the trajectory with h = h = 1/2 + i/^, which is obtained 
from the ground state trajectory with h = h = 1/2 + iuh by analytical continuation from real to 

^''We are grateful to J. Wosiek for making the Mathematica code used in (l4| available to us. 
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imaginary Uh- The same happens for the = solution at = 3 found there - its total spin is 
h = h = 1. As a consequence, the states found in the dVL paper do not belong to the quantum 
space of the SL{2, C) magnet and, therefore, are unphysical. 

We conclude that the criticism of our results in the dVL paper is groundless. At = 3 
the quantization procedure proposed in that paper does not reproduce correctly the part of the 
spectrum corresponding to Regs 7^ 0, whereas at A^ = 4 it generates spurious states, which do 
not belong to the quantum space of the model. 



A Appendix: Solution to the Baxter equation at = 2 



In this Appendix we summarize the properties of the eigenvalues of the Baxter operator Q{u, u) 
at N = 2. As was shown in the Section 3, Q{u,u) is equal to the integral of the function 
Q{z,z) defined in ( |3.32| ) over the two-dimensional plane. 

The function Qs{z\ h) entering (p.32|) is expressed in terms of the Legendre function on the 
second-kind, Eq. ( pj.31| ). Using the properties of the Legendre functions [Q, one finds the 
behaviour of Qs{z; h) around singular points 2; = and z = 1 as 



--lnz-^(l-/i)+^/'(l) + C(z) 



Qs{z-h) {l-z) 



2s-h-l 



r\i-h) 

_2T{2 - 2h) 



+ 0{l-z) 



(A.l) 
(A.2) 



Comparing ( |A.1| ) and ( |A.2| ) with the asymptotic behaviour of the functions Qn\z) and Qm{z) 
defined in ( p.7|) and (|3.16|) , one finds that Qs{z; h) can be decomposed over the fundamental set 
of solutions around z = as 



Qs{z\ h) = - 
Qs{z] l-h) 



^{l-h)-^{l)-- 



(0), 



and over the fundamental set of solutions around 2; = 1 as 

T\l-h) 



Qs{z; h) 



Q?{z), Qs{z-l-h) 



2V{2-2h) 



2T{2h) 



(A.3) 



(A.4) 



The functions Qn\z) and Qm'{z) are related to each other by the transition matrix Qn'{z) 



(0), 



'Ylfm^nmQni {z) defined in ( p.26| ). Comparison of the r.h.s. of (|A.3|) and ( |A.3|) yields 



n{h) 



A{h) A{l-h) 
S{h) 6{1 - h) 



(A.5) 



where A{h) = T{2h - l)/T^{h) and S{h) = -2 [^(/i) - ^(1) - |] A{h). Similar matrix in the 
antiholomorphic sector equals Q = Q(h). 

The general expression for the function Q{z, z) is given by ( |3.32D with the expansion coef- 
ficients c/i defined in ( p.33| ). Substituting ([A.3|) and ( |A.4| ) as well as analogous relations in the 
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antiholomorphic sector into (|3.32| ), one arrives at Eqs. ( p.9|) and (|3.20|) with the mixing matrices 
given by 



ai{h)-2 1 
1 



C(i) 



Ph 
Pi-h 



where ai{h) was defined in (|3.35|) and 



T\l-h)T{2h-l) 
' T^{h) V{2-2h)' 



(A.6) 



(A.7) 



It is straightforward to verify that the matrices C-"-* and C^^-* satisfy the quantization conditions 
To obtain the Q— block at = 2 one inserts ( g^STD into 



Q{u] h) 



T{2s - h) 



dz z"' 



'l-z 



Q 



-h 



l + Z 
l~Z 



Integration can be performed by replacing the Legendre function by its integral representation. 
In this way, one arrives at ( |4.20| ). The function Q{u,u) with required analytical properties is 
given by the bilinear combination of the holomorphic and antiholomorphic Q— blocks, Eq. ([4.21|) . 

As was explained in Section 4.4, to calculate the energy and quasimomentum it is convenient 
to introduce the blocks, Qo{u) and Qi{u), Eqs. ( |4.32|) and (|4.34|) , respectively. At A^ = 2 the 
block Qo{u) is fixed (up to an overall normalization) by the requirement to have simple poles at 
the points u = —i{s — m) with m > 0, ^ 



Qoiu) 



r(2s - 1 + h)T{2s - h) 
1 



dz z"' 



'l-z)2{-i)p_ 



l + z 
1 - z 



T^{2s] 



3-^2 



s — iu,2s — h,2s — 1 + h 
2s, 2s 



(A.9) 



The antiholomorphic block Qq{u) can be obtained from ( [A. 91 ) by replacing u —u, s —>■ s and 
h —>■ h. The block Qi{u) is defined according to (|4.34| ) as 



Qi{u) 



r(l -s + iu) 
r(s + iu) 



r(l -s + iu) 
r(s + iu) 



Qo{u] 1 - s) , 



(A.IO) 



where in the last relation we indicated explicitly the dependence of the block Qo on the spin s. 
The two sets of the blocks are linearly dependent 



1 



■Q{u;h) - 



1 



r(2s-l + M~^^^"^^ r(2s-/.)^^^''-^^ = -vrcot(vr/.)Qo(n 



T{l-2s + h) 
The inverse relation reads 



T{2-2s-h) 



7rcot(7r/i) Qi(m) . (A.ll) 



Q{u, h) = p(s, h) 



1 



Q{u, 1- h) = p{s, h) 



T{2-2s-h) 
1 



Qq{u) - 



r(2s - h) 



Qi{u) 



T{l-2s + h) 
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Qo{u) - 



T{2s-l + h) 



Qi{u) 



(A.12) 



where p{s, h) = ir^/ (2 sin(27rs) sm{nh)). Substituting (|A.12|) into ( [4.21| ) one finds that the result- 
ing expression matches ( |4.36|) . 

Using the properties of the 3F2— series [^], one can show that the block Qo{u), Eq. ( [A .91 ), 
satisfies the following relation 



sinfTr/i) „ , . „ 
Qo{-u; s) - ^— — -Qo{u; s) = -T 
sin(27rs) 



2s,l — 2s,l — s — iu,l — s + iu 
2s — h,2s — 1 + h, s — iu, s + iu 



Qo{u-l-s), (A.13) 



where we indicated explicitly the dependence on the spin s. Applying ( |A.1(J| ), ( |A.12| ) and ( |A.13| ) 
one can verify that the eigenvalue of the Baxter operator at = 2, Eq. ( |4.21|) , can be rewritten 
(up to an overall normalization) as 



Q{u,u) 



Qo{u) {Q,{-u*)r + (-1)"" Qo{-u) (Qo(r))* . (A.14) 



1 — s — in, 1 — s + m 
s — iu,s + iu 

Finally, two equivalent expressions for the eigenvalues of the Baxter Q— operator at A^ = 2, 
Eqs. ( A.14 ) and (|4.21|) , admit an elegant representation in terms of two-dimensional Feynman 
diagrams (see Figures 10a and b in Ref. |^). 



B Appendix: Properties of the Q— blocks 

In this Appendix we establish different useful relations between the blocks Q{u] h, q) and Q{u] h, q) 
defined in Eq. (|]T2|). 



Intertwining relations 

It is well-known that the SL{2, C) representations of the principal series of the spins (s, s) and 
(1 — s, 1 — s) are unitary equivalent. At the level of the eigenvalues of the Baxter Q— operator, 
this property leads to the following intertwining relation between the blocks 



Qs{u; h,q) 



r(l -s + iu) 
T(s + iu) 



N 



Qi-s{u] h,q) 



(B.l) 



where subscript indicates the corresponding value of the holomorphic spin. Indeed, one verifies 
that the both sides of this relation satisfy the Baxter equation (|2.7|) , have the same analytical 
properties ( [4.16 ) and asymptotic behaviour ( [4.17| ) at infinity. In a similar manner, using the 
identities s* = 1 — s, h* = 1 — h and qu = qli can show that 



{Qi-s{u*\ h, q))* = Qs{-u; I - h, -q) . 



(B.2) 



Combining together Eqs. ( p.l| ), ( |B.2|) and (|4.18| ), we obtain the relations between the holomorphic 
and antiholomorphic blocks 



Q{u;h,q) 
Q{u; h, q) 



r(l - s + iu) 
r(s -|- iu) 

r(l — s — iu) 
vis — iu) 



N 

{Qiu*;l-h,q)y 



N 



{Q{u*-l-h,q)y 



(B.3) 



Here, the ratio of F— functions compensates the difference in the analytical properties of two 
blocks and their asymptotic behaviour at infinity. 



49 



Wronskian relations 



The functions Q{u; h, q) and Q{u] 1 — h,q) satisfy the same Baxter equation (|2.?| ). This suggests 
to define their Wronskian as 



W{u) 



r{iu + s) 
T(iu — s) 



N 



Q{u + i] h, q) Q{u; 1- h,q) - Q{u; h, q) Q{u + h,q) 



(B.4) 



It follows from ( ^.71 ) that W{u) is periodic, W{u+i) = W{u). In addition, taking into account the 
properties of the blocks, Eqs. (f4.16|) and ( [4.17|) , one finds that W{u) is analytical in the half-plane 
Re{iu — s + 1) > and behaves there at large u as W{u) ~ u^. This implies that W{u) takes 
constant values for arbitrary u. Then, substituting ( |4.17| ) into (|B.4D , one finds W{u) = 1 — 2h . 
Similar consideration in the antiholomorphic sector leads to 



Q{u - i; h, q) Q{u; 1 - h,q) - Q{u; h, q) Q{u - i;l - h,q) = {1 - 2h) 



-lU 



N 



(B.5) 



For u = i{s + e) and u = —i{s + e) we find from the Wronskians 
$(e) and $(e), defined in ([4.26|) , satisfy the relations 

$(e)-$(e + A;) = 0(e^), 



T{—iu + s) 
and (|B.5[ ) that the functions 



<l>(e)-<l>(e + A;) = 0{e 



(B.6) 



with k being positive integer. Here, we used the fact that Q{i{s + e);h, q) and Q{i{s + e); 1 — /i, g) 
are finite for e — 0. 



Series representation 

We can obtain a series representation for the block Q{u;h,q) in the different regions on the 
complex M— plane, by replacing the function Qi{z) in ( |4.12| ) by its expressions in terms of the 
fundamental solutions, Eqs. ( [4.6| ) and ( [4.7| ), defined in ( |3.7| ) and ( |3.16| ). For z — > 1 the function 
Qi{z) is given by 



(B.7) 



n=0 



with t>o = 1 and the expansion coefficients Vn = Vn\—q) defined in ( |3.16| ) and ( |3.17| ). In this 
way, one gets from (^4.12 ) 



, . r(l-s + m)^ 



T{n-h + Ns) 



T{-h + Ns) T{n + 1 - h - s + iu + Ns) ' 



(B.8) 



For z —>■ the function Qi{z) is given by 



Q,{z)=z'-^J2(lnzrJ2 



(B.9) 



A;=0 



ra=0 



with Wn^ = C5_iMn and the expansion coefficients defined in ( |3.7| ) 

and (R.8|). This leads to 



^' " T(Ns - h) 5- 5^ Uu-s + n + 

^ ' n=0 k=0 ^ ' ' / 



(B.IO) 
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where the sum over n goes over the fc— th order poles located ai u = i{n + 1 — s). 

We notice that ( [B.8| ) reproduces correctly the asymptotic behaviour of Q{u; h, q) at large u, 
Eq. ( [4.17| ), and the position of its poles on the m— plane, Eq. ( |4.16| ), but not their order. The 
reason for this is that the series ( B.8 ) is convergent only for Re(l — s + iu) > 0. Indeed, the 
large-order behaviour of the expansion coefficients f„ in ( |B.7| ) is determined by the asymptotics 
of the function Qi{z) at ;z = 0, Eq. (|R9D 



1 

27ri 



dz 



Therefore, for Re(l 
^^ln^-^n/n2-^+^". 



|l-^|<e 

s + iu) 



< the series 



dz ■ 



In 



N-l 



In 



N-l 



n 



n 



.1 - zY+'^ 

diverges at large n as Xln '^n/'"'^"'^^™ 



B.ll) 



C Appendix: Contour integral representation 



Let us demonstrate that the two-dimensional integral ( |4.5|) can be decomposed into the sum 
of products of simple contour integrals, Eq. (|4.9|) . The derivation is based on the technique 
developed in |28| for calculation of the correlation functions in CFT. To simplify notations, we 



rewrite the integral over ( |4.5| ) as 



Q = J d^z ^ qn{z)Cnkqk{z) = j d^zq^{z) ■ C ■ q{ 



(C.l) 



Here, the functions qn{z) 
points located at Zi 



n,k=l 
-iu—l 



Qniz) and qkiz) 



Z 



Qk{z) have three isolated singular 



Zi = 0, 1 and oo. Around these points, they have a nontrivial monodromy 

(C.2) 



qn{z) ^ [Mi]nkqk{z) 



qn{z) ^ [Mi]nkqk{.z) 



with Mi and Mi being the corresponding monodromy matrices such that MqMiMoo = 1 and 
MqMiMoo = 1. Here, z and z encircle the singular points on the complex plane in the counter- 
clockwise and clockwise directions, respectively. For the integrand in ( |C.1| ) to be a single-valued 
function on the two-dimensional plane, the mixing matrix in Eq. (|C.1|) has to satisfy the condition 



MfCMi = C. 



0, l,oo). 



(C.3) 



Since q{z) and q{z) are analytical functions on the complex plane with two cuts running from 
the singular points and 1 to infinity, as shown in Figure one can apply the Stokes's theorem 



Q 



"'4 

oz 



q^z) . C 



dz q{z') 



2i 



- / dzq^{z)-C- 



as 



dz' q{z') 



(C.4) 



Here, z^nx is an arbitrary reference point and the ^'—integration goes along the contour that does 
not cross the cut. It is convenient to choose z^ux = and split the integral over (infinite) contour 
dT, into four integrals along the different edges of two cuts. Their contribution to ( p.4| ) can be 
calculated as follows 

Qlb'a'] = / dzq^{z)-C- I dz'q{z')+ I dz' q{z' 

J Qi _J 1 v 

Qh = ^jyzq^{z)-C- J^^ dz' qiz') + J\z' qiz') 



(C.5) 
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r 



00 



00 



b' 



9S 

Figure 10: The integration contour 9S in Eq. ( p.4|) . 

In these relations the integration goes along two different edges of the same cut. In spite of 
the fact that g(z) and q{z) are discontinuous across the cut, their bilinear combination remains 
continuous due to ( |C.3| ) 



h pz fb 



dzq^{z)C / dz'q{z')= / dz[Miq{z)YC / dz' Miq{z') = / dzq^{z)C / dz' q{z' 



This leads to a partial cancellation of terms in the sum Q[b'a'] + Q[ab]- Then, one calculates the 
contribution of the second cut 



Q[cd] = j dz q'^ (z) ■ C ■ dz' q{z' 

Q[d'c'] - 



dzq'{z)-C- dz'q{z') 
2* .Id' Jo 



(C.6) 



and finds that the same property leads to Q[cd\ + Qid'c'] = 0. Combining together (|C.5|) and (|C.6|) , 
we obtain the following expression for the two-dimensional integral ( P.4| ) 



Q 



1 

2i 



pb pb' pi 

/ dzq'^iz)- / dzq^{z) C- dz' q{z' 

Ja J a' Jo 

-t POO pi 

- / dz q^{z) ■ (1 - M^C ■ / dz' q{z') . 
2^ Ji Jo 



(C.7) 



Finally, one replaces the functions q{z) and q{z) by their actual expressions, Eq. ( |C.1| ), and 
arrives at ( [4.9| ). 

One can obtain another (though equivalent) representation for Q if one starts with 



Q= f (fz^ f dz'q^(z') ■ C ■ q(z) = -- [ dz f dz'q^iz') ■ C ■ qiz) 
is oz LJo J 2z Jq^ Jq 



instead of (|C.4| ). Repeating the same analysis one gets 



Q = - 



2i 



dz'q^iz') ■ C(l-Mi) 



dz q{z) 



(C.9) 
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Let us now take into account that q{z) is analytical inside S and, therefore, J^j^dzqlz) = 
0. As before, splitting the integration contour into four pieces and taking into account that 
Ml dz q{z) = j^, dz q{z) and dz q{z) = Mq J^f dz q{z) one obtains 

/oo poo 
dz q{z) + {1 - Ml) J dzq{z)=0. (C.IO) 

Here, two terms in the l.h.s. correspond to the contribution of two cuts. We conclude from ( P.10|) 
that 

pi poo 

-(l-Mo)y dz q{z) = {1 - MqMi) J dzq{z). (C.ll) 

Obviously, similar property holds for the function q{z) in the antiholomorphic sector. Recalling 
the definition of the function q{z), Eq. ( |C.1|) , one obtains 

1 rl roc pi 



dzqn{z)= rf^^-™-^g„(^), / dzqn{z)= dzz'^-'Qn{l/z). (C.12) 
^0 Jl Jo 

Eq. (|C.11| ) allows us to rewrite (|C.7|) and ( |C.9| ) in the form 

Q = jjyzq^{z) ■ (l - [l-M^y' - (1 - M^yy'c ■ j\-zq{z) , (C.13) 
in which the symmetry between z— and ^—sectors becomes manifest. 

D Appendix: Degenerate Q— blocks 

In our analysis of the quantization conditions performed in Section 4.3, we have tacitly assumed 
that the blocks Q{u\ h, q) are well-defined for arbitrary spins h = {1 + 71^) + and, in addition, 
Q{u] h, q) is finite at m = 2(5 + n — 1) with n > 1. As was already mentioned in Section 4.2, the 
first condition is not satisfied at Uh = 0. The second condition does not hold for (half)integer 
spins s since, by the definition, the block Q{u] h, q) has poles at u~ = —i{s — m) with m > 1 and 
for positive integer n and m, such that 2s — 1 = m — n, the point u = i{s + n—l) coincides with 
the pole M~ . In this Section, we will work out the quantization conditions for (half)integer spins 
h = {1 + nh)/'2, and s. It worth mentioning that one has to deal with these two cases calculating 



the ground state of the Schrodinger equation ( |1.2| ) at /i = 1/2 and s = 



To start with, let us examine the series representation (|B.7|) for the block Q{u\ h,q) ai h = 
(1 + nh)/2 + iufi in the limit zz/^ — and rih > 0. The expansion coefficients Vn entering ( [B.7| ) 
satisfy the iV— term recurrence relations, which lead to u„ ~ l/z//i as zz/^ ^ for n > 71^. The 
resulting expression for Q{u] h, q) can be written as 

Q{u- h, q) = ^^Q{u; l-h,q)+ Q{u; n^, q) + 0{ph) , (D.l) 

where Q{u] h, q) and Q{u; 1 — h,q) are finite for z//i and Anf^{q) = \im^i^^o[vn^{q)i^h]- Taking 
into account ([4.17|) , we find that the function Q{u;nh,q) defined in this way has the following 
asymptotic behaviour at large u 

Q{u;n,,q) ~ {zu)-''^+^'+^'-y' [1 + 0{l/u)] (D.2) 

+An,{q) (m)-^^+(i-"'')/2 \nu [1 + 0{l/u)] . 
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According to the definition (|D.1| ), the function Q{u] Uh, q) is a linear combination of two degen- 
erate blocks and, therefore, it satisfies the holomorphic Baxter equation (|2.7] ). Eqs. ( p . 1|) and 
( p.2| ) are valid only for rih > 0. At rih = 0, or equivalently h = 1/2, the function Q{u;0,q) is 
defined as 



Q{u] 0, q) = d^^Q{u] 1/2 + iuh, q) 



(D.3) 



It is straightforward to verify that for > the function Q{u; rih, q) satisfies the chiral Baxter 
equation ( ^^Tl) . 

Let us now insert ( p.l|) into the quantization conditions ( [4.28 ) and examine the limit ^ 0. 
It follows from (|;2|) that 



Q{i{s + e)]nh,q) 
' Q{i{s + e); (1 - nh)/2,q) 



+ 



:d.4) 



Then, comparing the coefficients in front of different powers of we find from (|4.28| ) the set of 
N — 1 quantization conditions on the charges ga, q^ 



Qn 



Im 



Q{i{s + e);nh,q) 



Q{t{s 



nh)/2,q) 



0, 



(D.5) 



e=0 



with n = 1, — 1 and h = {1 + nh)/2. We recall that the additional set of A^ quantization 
conditions follows from ( [4 .311) in the antiholomorphic sector. At /i = 1/2 the relation (p.5| ) leads 
to 

(9" (9 / ^ 

— \nQ\^i{s + e)]l/2 + iuh,q^ 







e=Uh=0 



(D.6) 



Let us now consider the quantization conditions for (half) integer spins s. In distinction with 
the previous case, the blocks Qs{u] h, q) as well as the eigenvalues of the Baxter operator Q{u, u) 
remain finite in the limit s = (1 + ns)/2 and possess a correct analytical properties as functions 
of u and u. Nevertheless, the quantization conditions have to be modified because the two sets 
of points ( [4.25|) overlap. Repeating the analysis of Section 4.4, we find that the function Q{u,u) 
has correct analytical properties for (half)integer s provided that the following conditions are 
satisfied 

Q{i{l+\ns\)/2 + e;h,q) 



arg 



©.,5 + Oie 



2N\ 



(D.7) 



[Q{t{l + \ns\)/2 + e;l-h,q)\ 

Here, in distinction with ( [4.30|) , the Q— blocks are calculated in the vicinity of the point u = + 
\ns\)/2 that belongs to the both sets ( [4.25| ) simultaneously. In addition, the small e— expansion 
in the r.h.s. starts with the terms ~ e^^. 
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